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ERROR  PERFORMANCE  OF  DIRECT  SEQUENCE  SPREAD 
SPECTRUM  SYSTEMS  ON  NONSELECTIVE 
GENERALIZED  RAYLEIGH  FADING  CHANNELS 

Robert  Charles  Hanlon,  Ph.D.,  Major,  USAF 
Department  of  Electrical  Engineering 
University  of  Illinois  at  Urbana-Champaign ,  1980 

Abstract 

The  error  performance  of  direct -sequence  spread  spectrum  multiple-access 
communications  systems  on  nonselective  fading  channels  is  investigated.  The 
channels  considered  are  those  for  which  the  fading  intensity  is  slow, 
nonselective,  and  exhibit  a  generalized  Rayleigh  distribution.  „The  generalized 
Rayleigh  distribution  encompasses  the  important  special  cases  of  Rician 
(specular  plus  Rayleigh)  fading,  Nakagami  m-distributed  fading,  and  Rayleigh 
fading.  Results  are  obtained  for  the  special  cases  of  a  single  cochannel 
interferer  without  additive  noise  and  multiple  independent  cochannel 
interferers  with  additive  white  Gaussian  noise.  The  results  are  developed  as  a 
series  expansion  in  the  moments  of  the  cross-correlation  between  the  spreading 
codes  (single  interferer)  or  as  a  series  expansion  in  the  powers  of  the  second 
moment  of  the  code  cross-correlation  (multiple  interferers).  To  develop  simple 
estimates  of  the  probability  of  error,  empirical  estimates  of  the  low-order 
cross-correlation  moments  of  the  spreading  codes  (taken  to  be  Gold  codes)  are 
developed  from  a  generalized  gamma  density  approximation  based  on  the 
cross-correlation  moments  for  both  Gold  codes  and  random  codes.  Numerical 


results  for  the  probability  of  error  are  obtained  with  these  estimates  and 
compared  with  the  results  based  on  exact  moments  developed  from  the  code 
cross-correlation  probability  density.  (103  pages). 
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Chapter  1 
INTRODUCTION 

1.1.  Spread  Spectrum  Communications  S,ys,tSB3 

A  spread  spectrum  communications  system  is  characterized  by  the  use  of  much 
larger  radio  frequency  bandwidths  than  required  by  the  bandwidth  of  the 
baseband  information  signal  [Dixon,  1976a],  [Glazer,  1973].  By  effective  use 
of  this  additional  bandwidth,  spread  spectrum  systems  develop  several  inherent 
advantages  over  systems  which  use  conventional,  relatively  narrowband 
techniques.  These  advantages  include  [Cahn,  1973],  [Dixon,  1976a],  [Glazer, 
1973]:  (1)  enhanced  resistance  to  interference,  (2)  capability  for 
multiple-access  within  a  common  frequency  spectrum,  (3)  ability  to  operate 
effectively  within  a  frequency  spectrum  simultaneously  being  occupied  by 
conventionally  modulated  signals,  (4)  low  detectibility  by  unintended 
receivers,  and  (5)  message  privacy.  These  benefits  are  interrelated  and  arise 
due  to  the  methods  employed  to  spread  the  bandwidth  prior  to  transmission  and 
compress  (despread)  the  bandwidth  after  reception  by  the  intended  recipient. 
However,  these  techniques  also  give  rise  to  two  principal  disadvantages: 
relatively  complicated  transmitting  and  receiving  equipment,  and  conflict  with 
current  frequency  allocation  and  assignment  procedures. 

There  are  several  methods  which  are  commonly  used  to  spread  the  information 
signal  bandwidth.  These  include  direct  sequence  code  modulation,  frequency 
hopping,  time  hopping,  swept  frequency  ("chirp")  modulation,  and  hybrid 
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combinations  of  these  methods.  The  basic  underlying  principle  of  spread 
spectrum  systems  is  that  the  intended  recipient  knows  both  the  technique 
employed  and  its  specific  implementation.  The  intended  recipient  is  then  able 
to  use  this  knowledge  to  invert  the  spreading  process,  thereby  despreading  the 
received  signal  to  reproduce  the  narrowband  information  signal.  To  clarify 
these  concepts,  direct  sequence  modulation  is  examined  in  detail.  Further,  the 
scope  of  the  remainder  of  this  paper  is  restricted  to  this  form  of  spread 
spectrum  modulation. 

1.2.  Direct  Seausnse  Code  Modulation 

In  a  direct  sequence  (DS)  code  modulation  system,  the  information  signal 

(digital  or  analog)  is  multiplied  by  a  spreading  code  sequence  whose  effective 

bandwidth  is  much  greater  than  that  of  the  information  signal.  This 

multiplication  in  the  time  domain  is  equivalent  to  convolution  in  the  frequency 

domain.  If  the  baseband  information  signal  is  nominally  bandlimited  to  the 

interval  (-Bb,  Bb)  and  the  spreading  code  sequence  is  nominally  bandlimited  to 

the  interval  (-Bm,  Bm) ,  then  the  product  signal  in  general  occupies  the 

frequency  interval  (-B.-B  ,  B.+  B  ).  Since  in  spread  spectrum  systems,  B  is 

d  m  d  m  m 

much  greater  than  B^,  then  the  product  signal  bandwidth  is  approximately  that 
of  the  spreading  code  sequence.  For  high  rate  spreading  codes,  the  transmitted 
power  spectral  density  can  be  made  quite  low.  As  a  result  of  the  low  power 
spectral  density,  interference  effects  on  conventionally  modulated  systems 
operating  simultaneously  over  some  small  portion  of  the  spread  spectrum 
bandwidth  can  be  minimal,  and  the  detectibility  of  the  spread  spectrum 
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transmission  by  conventional  receivers  is  significantly  reduced.  Low 
detectibility  enhances  message  privacy  and  traffic  flow  security  (precludes 
knowledge  of  message  traffic  volume)  and  compounds  problems  in  radio  direction 
finding. 

For  the  purposes  of  this  discussion,  it  is  assumed  that  the  information 

signal  is  a  binary  digital  signal,  b(t),  and  the  spreading  code  sequence 

(sometimes  referred  to  as  a  key  sequence  or  signature  sequence)  is  also  a 

binary  digital  signal,  a(t),  where  both  signals  can  take  the  value  of  either 

plus  or  minus  one  during  their  respective  bit  intervals.  The  bit  interval  of 

the  information  signal  is  Td  and  the  bit  interval  for  the  spread ins  code 

sequence  is  Tc,  where  Td  is  much  greater  than  Tc  to  provide  the  required 

greater  bandwidth  for  the  spreading  code.  A  bit  of  the  spreading  code  is 

referred  to  as  a  chip  to  distinguish  it  from  an  information  signal  bit, 

consequently,  T  is  referred  to  a3  the  chip  interval  and  T.  as  the  bit 

c  d 

interval.  The  spreading  code  is  repetitive  with  period  LTq,  where  L  is  the 
number  of  chips  per  period,  and  is  referred  to  as  the  code  length.  It  is 
assumed  that  the  spreading  code  is  synchronized  to  the  information  signal  such 
that  one  period  of  the  spreading  code  occurs  during  each  bit,  i.e., 


(1.1) 
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The  product  signal  is  Ca(t)b(t)].  If  the  product  signal  is  multiplied  by 

p  2 

a(t),  we  have  [a  (t)b(t)].  Since  a(t)  is  either  plus  or  minus  one,  a  (t)  is 

2 

identically  one,  and  [a  (t)b(t)]  is  equal  to  b(t),  which  is  the  original 
information  signal.  Therefore,  the  despreading  operation  is  identical  to  the 
spreading  operation,  that  is,  multiplication  of  a  signal  by  the  spreading  code. 
Clearly,  the  receiver  requires  prior  knowledge  of  the  spreading  code.  Hence, 
the  spreading  code  must  be  deterministic.  However,  if  the  code  length  is 
fairly  long,  the  spreading  code  appears  to  be  random  and  is  often  referred  to 
as  a  pseudorandom  code  sequence.  The  receiver  must  not  only  know  the  code 
sequence,  but  must  also  generate  the  sequence  in-phase  with  the  sequence  which 
has  been  applied  to  the  information  signal.  In  a  system  employing  radio 
propagation,  a  pure  synchronous  system  is  generally  not  possible.  This  is  the 
reason  for  using  a  periodic  sequence.  The  periodic  nature  of  the  code  sequence 
enables  the  receiver  to  detect  the  phase  and  to  synchronize  to  thi3  phase. 

/ 

If  signals  other  than  the  intended  spread  spectrum  signal  are  received, 
multiplication  by  the  spreading  code  in  the  receiver  spreads  each  of  these 
undesired  signals  to  at  least  the  bandwidth  of  the  spreading  code.  Narrowband 
filtering  after  the  multiplication  passes  the  despread  intended  signal  while 
rejecting  most  of  the  power  of  unintended  signals.  This  can  be  seen  to  be  a 
correlation  receiver  where  the  correlation  is  against  the  spreading  code  of  the 
intended  signal.  For  conventional  narrowband  signals,  their  power  spectral 
densities  are  reduced  by  a  factor  of  approximately  (Bb/Bm)«1  by  the  receiver's 
despreading  process.  Even  if  the  interference  has  the  same  bandwidth  as  the 
spread  spectrum  signal,  its  power  spectral  density  would  be  reduced  by 
approximately  (Bn/2Bm)=1/2,  that  is,  3  dB.  Further,  a  signal  whose  bandwidth 


was  on  the  order  of  the  spread  spectrum  bandwidth  would  generally  have  an 
inherently  low  power  spectral  density.  To  help  minimize  noise  and  interference 
effects,  a  wideband  filter  can  be  used  on  the  received  signal  prior  to 
despreading  to  eliminate  undesired  components  outside  the  spread  spectrum 
bandwidth.  It  has  also  been  shown  that  notch  filters  can  be  used  to  eliminate 
narrowband  interference  wihout  undue  distortion  of  the  spread  spectrum  signal 
[Sussman  and  Ferrari,  1976]. 


1.3.  Direct  Sequence  Multiple-Access  Communications 

As  indicated,  the  despreading  process  is  fairly  effective  against  undesired 
signals.  This  feature  can  be  used  to  obtain  multiple-access  capability  by 
assigning  each  of  the  various  system  users  a  unique  spreading  code.  Provided 
that  the  codes  used  do  not  exhibit  a  high  degree  of  cross-correlation,  then  the 
intended  receiver  is  able  to  effectively  reject  the  signals  from  the  other 
users.  Use  of  codes  to  implement  multiple-access  communications  is  referred  to 
as  code  division  multiplexing  and  offers  an  attractive  alternative  to 
conventional  time  and  frequency  division  multiplexing.  The  presence  of  the 
other  users  in  code  division  multiplexing  results  in  an  inherent  cochannel 
interference.  To  help  minimize  the  cross-correlation  between  the  assigned 
spreading  codes  and  thereby  reduce  the  cochannel  interference,  special  codes 
such  as  Gold  codes  are  generally  used  which  are  known  to  have  good 
cross-correlation  properties  [Di>...,  1976a],  [Gold,  1964],  [Pursley,  1976].  It 
should  be  noted  that  when  Gold  codes  or  other  linear  codes  are  used,  they  do 
not  provide  message  security  in  the  sense  of  cryptographic  codes  [Dixon, 


1976a],  [Scholtz,  1977],  [Ristenbatt  and  Daws,  1977],  but  do  provide  a  degree 
of  message  privacy. 

Figure  1  shows  a  basic  direct  sequence  spread  spectrum  multiple-access 
(DS/SSMA)  communications  system  [Gardner  and  Orr,  1979],  [Orr,  1977].  The 
receiver  is  assumed  to  be  frequency,  phase  and  code  synchronized  to  the 
intended  signal  (channel  1).  Although  acquisition  and  maintenance  of 
synchronization  is  a  major  aspect  of  spread  spectrum  communications  systems,  it 
is  beyond  the  scope  of  this  paper  and  is  assumed  to  be  perfect  in  all 
subsequent  analyses.  Although  the  receiver  is  assumed  to  be  synchronized  to 
channel  1 ,  the  transmitters  are  operating  asynchronously  relative  to  each 
other. 

1.4.  Channel  Model 

In  analyzing  the  performance  of  a  communications  system,  the  channel  over 
which  the  signals  propagate  must  be  modeled  to  account  for  all  significant 
sources  of  signal  degradation.  Along  with  the  presence  of  additive  random 
noise,  it  is  relatively  common  in  systems  utilizing  radio  propagation  to 
encounter  random  multiplicative  signal  distortion,  that  is,  continuously  time 
varying  received  signal  amplitudes.  This  random  multiplicative  signal 
distortion  is  generally  referred  to  as  nonselective  fading.  Fading  generally 
arises  out  of  conditions  which  give  rise  to  multipath  propagation  (e.g., 
ionospheric  reflection  or  tropospheric  scatter  propagation)  where  the  multiple 
paths  are  of  varying  lengths  and  nonstationary  in  nature  (e.g.,  dynamic 
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inhomogeneities  in  scattering  volume).  Nonselective  fading  implies  that  the 
fading  effects  do  not  change  appreciably  with  frequency  over  the  bandwidth  of 
interest.  In  addition  to  assuming  that  the  fading  is  nonselective,  it  is  also 
assumed  that  the  fading  is  slow  ,  that  is,  the  fading  intensity  does  not  change 
appreciably  over  time  periods  on  the  order  of  the  data  bit  interval.  For  the 
purposes  of  this  analysis,  the  additive  random  noise  is  assumed  to  be 
independent  zero-mean  white  Gaussian  noise  (uniform  two-sided  power  spectral 
density  of  nQ/2  with  normally  distributed  amplitude).  Figure  2  depicts  the 
channel  model  [Gardner  and  Orr,  1979]. 

The  fading  intensities,  f^t),  are  assumed  to  be  generalized  Rayleigh 
variates  (Appendix  B).  The  generalized  Rayleigh  distribution  encompasses  a 
wide  class  of  distributions  which  are  useful  in  the  study  of  communications 
over  fading  channels.  In  particular,  the  generalized  Rayleigh  distribution  has 
as  special  cases:  the  Rician  (specular  plus  Rayleigh),  Nakagami 
(m-distributed) ,  and  Rayleigh  distributions.  Each  of  these  distributions  has 
been  extensively  utilized  for  modeling  fading  channels. 


1.5.  System  Equations 

For  the  spread  spectrum  multiple-access  communications  system  depicted  in 
Fig.  1,  the  transmitted  signal  on  the  i-th  channel  can  be  represented  as  an 
amplitude  modulated  carrier,  that  is,  [Gardner  and  Orr,  1979] 


a.(t)  b.(t)  cos  u  t 
i  i  c 


Si(t)  =  (2PTi) 


1/2 


(1.2) 


propagation  delay,  0^  is  the  random  carrier 
)  is  l he  additive  white  Gaussian  noise. 
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where  PTi  is  the  transmitted  power,  a^(t)  is  the  spreading  code,  b^t)  is  the 
information  signal,  and  coc  is  the  common  carrier  frequency.  The  channel  model 
as  illustrated  in  Fig.  2  introduces  nonselective  random  fading  and  random  path 
delays  on  each  signal  and  random  additive  noise.  The  total  received  signal  for 
the  case  of  M  users  is  then  [Gardner  and  Orr,  19791 


*(t)  =  n(t)  +  1  (2P . ) 1/2  f  (t-r  )  a  (t-r  ) 

i=  1  1  1  11  1 


b.  (t-r. )  cos  (<o  t  +  0, ) 
11  cl 


(1.3) 


where  the  P^  are  the  effective  received  powers  in  the  absence  of  fading,  the 
fi  are  the  random  fading  intensities,  the  are  random  time  delays,  and  the 
9^  are  random  phase  angles.  The  f^,  t\,  and  0.^  are  random  processes;  however, 
the  time  dependence  is  not  explicitly  shown  to  simplify  the  notation.  The  t 
are  assumed  to  be  independent  and  identically  distributed  (i.i.d.)  with 
uniform  distributions  on  the  interval  (0,  T^).  The  0^^  are  assumed  to  be 
i.i.d.  with  uniform  distributions  on  the  interval  (0,  2ir).  The  assumption  that 
the  9,^  are  uniform  implies  the  transmitters  are  operating  asynchronously  (as 
previously  mentioned)  and  that  the  transmitters'  relative  phase  drifts 
overshadow  their  phase  variations  due  to  fading  effects. 


For  coherent  reception  and  under  the  stated  assumption  that  the  intended 
receiver  is  perfectly  synchronized  (frequency,  phase,  and  spreading  code)  to 
the  first  signal,  detection  is  accomplished  by  multiplying  the  received  signal 
by  [a^t-T^  cos(&>ct+©1 )]  and  integrating  over  a  data  bit  period,  that  is, 
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correlating  the  received  signal  with  the  spreading  code  of  the  intended  signal. 
Since  coherent  detection  is  assumed,  we  can  set  T1  equal  to  zero,  ®  equal  to 
zero,  and  restrict  our  interest  to  the  data  bit  received  during  the  interval 
(0,  T^)  without  loss  of  generality.  The  carrier  frequency  is  much  greater  than 
T”1 ,  hence  the  contribution  from  the  double  frequency  components  is  negligible 
and  the  output  of  the  correlation  integrator  is 


R*b.  S  +  I+  H 


(1.4) 


where  the  intended  signal  component  is 


T 

td 


s  r  (p./2)uz  /  at  f.(t)  , 
1  o  ' 


(1.5) 


the  cochannel  interference  component  is 


T 

I  *  ?  (P  /2)1/2  /  dt  a  (t)a1(t-Ti)b1(t-Ti)fi(t-T1)cos  S*  , 
i=2  0 


(1.6) 


and  the  additive  noise  component  is 


N  =  J  dt  a.(t)  n(t)  cos  “  t 
0  1  c 


(1.7) 


Under  the  assumption  of  slow  fading,  f 1 ,  f^,  and  0^  are  relatively  constant 


over  the  integration  interval.  We  can  then  write 
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S  =  (P,T 2/2)1/2  f, 
i  a  i 


(1.8) 


and 


M 

■  l 

1=2 


(P^/2)1'2 


*11 


f^  COS  0^^ 


(1.9) 


where  c  is  the  normalized  cross-correlation  between  codes  a,  and  a.  in  the 
n1i  1  1 

presence  of  data,  that  is, 


c 


n 


1i 


*d 

/ 

0 


dt  a1(t)ai(t-Ti)bi(t-Ti)  . 


(1.10) 


For  equally  likely  data  bits,  the  detection  threshold  is  zero  due  to  symmetry, 
and  the  probability  of  error  is 


Pr{E)  =  g  [Pr{S+I-*-N<0|b1=1)  PH-S+I+JDOjb^-l}]  .  (1.11) 

By  letting  (Z  =  I+N)  and  by  noting  that  the  probability  density  of  Z  is 

symmetric,  we  obtain  for  the  average  probability  of  error 

oo  CU 

Pr(E}  =  /  dco  /  dv  p  (v,«)  s  —  Pr{  Iz  I  >  S)  (1.12) 

0  0  sz  * 


where  p  is  the  joint  probability  density  of  the  random  processes  S  and  Z. 

3Z 
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The  expression  for  the  probability  of  error  given  by  Eq.(1.12)  appears 

quite  simple;  however,  determination  of  the  joint  probability  density  p  and 

sz 

subsequent  evaluation  of  the  integrals  is  in  general  quite  complicated  due  to 
the  numerous  random  processes  involved.  Three  general  approaches  have  been 
used  to  evaluate  or  approximate  the  system  performance.  The  first  involves 
determination  of  the  equivalent  signal- to-noise  ratio  (SNR)  resulting  from  the 
cochannel  interference  and  additive  noise,  that  is,  [Gardner  and  Orr,  1979], 
[Pursley,  1977] 

SNR  =  <S2>  /  Var(Z)  .  (1.13) 

The  equivalent  SNR  can  then  be  taken  as  an  indirect  measure  of  the  system 
performance  (bit  error  rate)  or  an  estimate  of  the  probability  of  error  can  be 
developed  by  using  the  equivalent  SNR  in  the  well-known  error  expression  for 
signals  in  Gaussian  noise 


Pr{E}  =  ^  erfc[(SNR/2)1/2]  (1.14) 

where  erfc(x)  is  the  complementary  error  function.  The  second  approach 
involves  bounding  the  error  probability  by  utilizing  the  theory  of  moment 
spaces  [Borth,  1979],  [Yao,  1976  and  1977].  The  third  approach,  and  that 
addressed  in  the  remainder  of  this  paper,  is  to  develop  a  series  expansion  for 
the  error  rate  in  the  cross-correlation  moments  of  the  spreading  codes  [Gardner 
and  Orr,  1979],  [Hanlon  and  Gardner,  1979  and  1980].  Provided  that  the  series 


expansion  converges  rapidly,  the  first  few  terms  provide  an  accurate  estimate 


of  the  average  error  rate 


1.6.  Outline  tils.  Thesis 

As  indicated  in  the  preceding  section,  the  approach  used  in  this  paper  to 
evaluate  the  error  performance  is  the  use  of  series  expansions  of  the  error 
rate  in  the  cross-correlation  moments  of  the  spreading  codes.  To  utilize  this 
approach,  we  require  access  to  the  code  cross-correlation  moments.  This 
problem  is  addressed  first. 

As  seen  later,  in  general,  evaluation  of  the  cross-correlation  moments 
requires  knowledge  of  the  cross-correlation  probability  density.  Obtaining  the 
probability  density  can  be  very  difficult,  and  to  circumvent  this  problem, 
accurate  estimates  of  the  moments  are  desired.  In  Chapter  2,  an  algorithm  is 
developed  which  can  be  used  to  evaluate  the  probability  density  from  the  joint 
probability  mass  function  for  adjacent  cross-correlation  values  for  integral 
offset.  This  algorithm  is  used  to  evaluate  the  moments  for  the 
cross-correlation  for  a  set  of  Gold  codes  of  length  31  and  a  subset  of  Gold 
codes  of  length  127.  Gold  codes  are  used  throughout  the  analysis  as 
representative  spreading  codes  for  illustration.  These  moments  and  the  moments 
of  random  codes  of  the  same  length  (Appendix  D)  are  used  to  develop  an 
empirical  estimate  for  the  low-order  cross-correlation  moments  for  Gold  codes 
based  on  modeling  the  cross-correlation  distribution  by  a  generalized  gamma 
distribution.  The  empirical  estimates  developed  from  the  gamma  density  model 
are  used  in  the  estimates  of  the  probability  of  error  developed  in  Chapter  3. 


The  error  performance  analysis  in  Chapter  3  addresses  two  special  cases. 
First,  the  case  of  a  single  cochannel  interferer  in  the  absence  of  additive 


noise.  This  case  represents  a  performance  bound  for  multiple-access 
communications,  that  is,  it  represents  an  irreducible  error  rate  in  the 
presence  of  cochannel  interference  and  fading  for  arbitrarily  large  signal- 
to-Gaussian-noise  ratios.  The  second  case  addressed  is  for  multiple 
independent  interferers  with  additive  Gaussian  noise.  This  case  also  gives 
rise  to  an  irreducible  error  rate  for  arbitrarily  large  signal-to-Gaussian- 
noise  ratios.  As  part  of  the  analysis,  simple  single-term  estimates  for  these 
irreducible  error  rates  are  developed  and  evaluated. 

In  both  cases  mentioned  above,  the  analytical  results  are  based  on  the 
developments  given  in  Appendix  B  for  generalized  Rayleigh  variates.  In 
particular,  the  probability  of  error  is  expressible  in  terms  of  the  expressions 
for  the  probability  that  one  generalized  Rayleigh  variate  exceeds  another  for 
either  independent  or  correlated  variates. 

Chapter  4  summarizes,  provides  conclusions  developed  from  the  results 
presented  earlier,  and  gives  recommendations  for  further  study  in  this  area. 

Throughout  the  paper,  reference  is  made  to  mathematical  functions  and 
formulas  which  are  listed  in  Appendix  A. 


A 
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Chapter  2 

CROSS-CORRELATION  BETWEEN  SPREADING  CODES 

2.1.  Cross-Correlation  In  the  Presence  Data 

In  a  direct  sequence  spread  spectrum  communications  system,  multiple-access 
capability  within  a  common  frequency  spectrum  is  obtained  by  assigning  each 
system  user  a  unique  periodic  spreading  code.  The  presence  of  multiple  users 
within  a  shared  frequency  spectrum  results  in  an  inherent  cochannel 
interference.  However,  if  the  cross-correlation  between  the  spreading  codes  is 
kept  low,  then  the  cochannel  interference  effects  are  relatively  minimized 
(subject  to  constraints  such  as  number  and  proximity  of  other  users)  and  the 
likelihood  of  reliable  communications  is  enhanced. 

For  direct  sequence  multiple-access  systems,  the  output  of  the  correlation 
receiver,  which  is  assumed  to  be  synchronized  to  the  intended  signal,  contains 
cochannel  interference  terms  of  the  form  [Gardner  and  Orr,  1979],  [Pursley, 
1977] 


/  dt  a  ,(t-  )  a.  (t-T.-®. )  b.(t-x. ) 

q  j  J  ^  lii  i 


(2.1) 


where  is  the  spreading  code  on  the  intended  signal,  a^  and  b.^  are  the 

spreading  code  and  data  bit  respectively  on  the  i-th  interfering  user's  signal, 
and  (Jk  set  the  phase  of  the  spreading  codes  relative  to  their  data  bits,  L 
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is  the  length  of  the  spreading  codes  (number  of  chips  in  one  period  of  the 
code),  Td  is  the  data  bit  interval  (assumed  to  be  identical  for  all  users),  and 
is  the  relative  timing  offset  in  the  receiver  between  the  intended  user's 
and  interfering  user's  signals.  This  timing  offset  is  due  to  asynchronism 
between  transmitters  as  well  as  differing  propagation  delays  over  the  channel. 
The  are  assumed  to  be  uniformly  distributed  on  the  interval  (O.T^).  Both 

the  spreading  codes  and  data  signal  can  assume  values  of  +1  or  -1.  The 
spreading  codes  have  a  period  of  L  chips  where  each  chip  is  of  duration  Tc< 
This  period  is  assumed  to  correspond  to  one  data  bit  interval. 

The  integral  in  Eq.(2.1)  is  similar  to  the  cross-correlation  between  the 
two  spreading  codes  except  for  the  presence  of  the  interfering  signal's  data 
bit.  Consequently,  it  is  referred  to  as  the  cross-correlation  in  the  presence 
of  data  and  designated  c.  The  factor  b^  considerably  complicates  analysis  of 
such  forms.  Although  special  sets  of  spreading  codes  which  have  good 
cross-correlation  properties  (typically  meaning  bounded  by  some  value  much  less 
than  L)  are  generally  used,  the  presence  of  the  data  bit  significantly  alters 
these  properties  [Massey  and  Uhran,  1975].  For  instance,  the  magnitude  of 
cros3-correlation  values  for  Gold  codes  [Gold,  1964]  of  length  31  never  exceeds 
9.  However,  with  data  present,  the  magnitude  of  the  integral  in  Eq.(2.1)  for 
the  same  codes  can  be  a3  large  as  21  [Orr,  1977],  [Massey  and  Uhran,  1975]. 

In  evaluating  the  multiple-access  system  performance,  two  limiting  cases 
are  of  particular  interest:  the  case  of  a  single  interferer  and  the  case  where 
there  are  several  simultaneous  inter ferers.  When  referring  to  a  single 
interferer,  this  does  not  imply  that  there  is  only  one  other  system  user,  but 
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rather  that  only  one  of  the  other  system  users  is  simultaneously  transmitting 
(this  case  would  typically  occur  if  the  various  users  each  had  a  low 
transmission  duty  cycle).  Likewise,  in  the  case  of  several  interferers,  3ay 
M-1 ,  there  are  in  general  more  system  users  than  M.  However,  there  are  only  M 
system  users  (M-1  interferers)  simultaneously  transmitting.  In  the  case  of 
several  interferers,  the  composite  interference  is  a  sum  of  integrals  of  the 
form  given  by  Eq.(2.1).  Since  both  the  data  sequences  and  timing  offsets  are 
independent,  the  interference  in  the  multiple  interferer  case  tends  toward  a 
Gaussian  random  process.  In  either  case,  the  probability  of  error  is  a 
function  of  c.  To  exactly  evaluate  the  error  rate,  the  probability  density  of 
c  would  be  needed.  However,  as  shown  later,  determination  of  the  probability 
density  can  be  quite  complicated.  For  a  single  interferer,  the  problem  of 
determining  the  probability  density  can  sometimes  be  circumvented  by  expanding 
the  conditional  probability  of  error  given  c  in  a  series  involving  the  powers 
of  (c/L).  An  estimate  of  the  error  rate  can  then  be  obtained  by  truncating  the 
series  and  taking  the  expectation.  Thi3  approach  has  been  shown  to  be  quite 
useful  for  estimating  error  rates  of  direct  sequence  systems  operating  over 
nonselective  slow  fading  channels  [Gardner  and  Orr,  1979],  [Hanlon  and  Gardner, 
1979  and  1980].  For  the  case  of  multiple  interferers,  if  we  approximate  the 
interference  by  a  Gaussian  random  process,  we  can  often  restrict  our  interest 
to  determination  of  the  second  moment  of  the  cross-correlation.  Therefore,  in 
many  cases  the  probability  of  error  can  be  accurately  estimated  by  determining 
the  first  few  moments  of  c  or,  more  practically,  an  estimate  of  these  moments. 
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For  a  given  multiple-access  communications  system,  the  code  that  is 
designated  a..  depends  on  which  user  is  taken  to  be  the  intended  user  for  the 
analysis.  Likewise,  depends  on  the  signal  whose  contribution  to  the 
cochannel  interference  is  to  be  determined.  In  general,  it  could  be  any  of  the 
other  users.  From  a  systems  viewpoint,  we  are  concerned  with  the  performance 
of  the  system  over  all  possible  combinations  of  users.  Further,  the  codes 
assigned  to  any  or  all  users  might  be  occasionally  changed  due  to  operational 
considerations  (e.g.,  use  of  dynamic  code  assignments  to  enhance  the  long-term 
privacy  or  antijam  characteristics).  Consequently,  exact  analysis  for  a 
particular  pair  of  codes  is  of  limited  utility,  and  we  are  concerned  primarily 
with  the  expected  performance  (average  error  rate)  over  the  entire  set  of 
permitted  codes.  The  analysis  then  addresses  the  average  over  all  possible 
users  (i.e.,  aj  and  a.^  are  assumed  to  be  equally  likely  to  be  any  of  the 
possible  code  sequences,  where  i  is  not  equal  to  j).  This,  coupled  with  the 
presence  of  the  randomizing  effects  due  to  the  interfering  signal's  data 
sequence  and  the  relative  timing  offsets  (t^),  leads  to  the  conjecture  that  the 
low-order  moments  of  the  average  cross-correlation  in  the  presence  of  data  is 
close  to  those  for  random  code  sequences  of  the  same  length  (L  chips  per  data 
bit).  This  conjecture  is  heuristically  tested  by  example. 

In  Section  2.2  an  algorithm  for  calculating  the  probability  density  for  the 
code  cross-correlation  in  the  presence  of  data  is  developed.  Although  our 
immediate  interest  is  in  determining  the  average  over  all  possible  users 
(entire  code  set),  the  algorithm  can  be  directly  applied  to  any  subset  of  the 
codes  if,  in  a  specific  application,  a  priori  considerations  warrant  such  a 
restriction.  This  probability  density  can  then  be  used  to  obtain  exact  moments 


for  whatever  code  set  or  subset  is  under  consideration.  The  probability 


density  is  used  to  determine  the  moments  of  the  code  cross-correlation  in  the 
presence  of  data  for  a  set  of  Gold  codes  of  length  31  and  for  a  subset  of  Gold 
codes  of  length  127. 

In  Section  2.3  two  approaches  for  estimating  the  low-order  Gold  code 
moments  are  addressed.  The  first  method  directly  estimates  the  low-order 
moments  by  the  corresponding  moments  for  random  code  sequences  of  length  L 
(Appendix  D).  The  second  method  empirically  refines  these  estimates  by  using 
the  moments  of  a  generalized  gamma  density  for  which  the  parameters  are  chosen 
to  provide  close  agreement  with  the  Gold  code  moments. 


In  Eq.(2.1),  the  product  a^a^b^  is  not  periodic  due  to  the  presence  of  the 
factor  b^.  However,  since  the  correlator  is  restricted  to  the  interval 
(0,Td),  the  integrand  can  be  considered  to  be  a  sample  of  the  product  of 
periodic  codes  defined  respectively  by  a^  and  the  product  a^^  on  the 
integration  interval.  The  cross-correlation  between  a  pair  of  periodic 
synchronous  codes  is  given  by  the  number  of  agreements  minus  the  number  of 
disagreements  when  the  codes  are  compared  chip-by-chip  over  one  period  of  the 
codes  [Gold,  1964],  [Dixon,  1976a].  For  the  codes  to  be  synchronous,  we 
require  the  codes  to  have  the  same  period  and  the  same  chip  intervals  and  for 
the  chip  intervals  to  coincide  in  time.  For  a  pair  of  codes  of  length  L 
(period  of  L  chips),  we  have 
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c  =  2A  -  L  (2.2) 

where  c  is  the  cross-correlation  value  and  A  is  the  number  of  agreements. 
Hence,  c  is  even  for  L  even  and  odd  for  L  odd  (c  mod2  =  L  mod2)  and  is 
restricted  to  the  closed  interval  [-L,  L]. 

For  and  the  product  a^^  to  be  synchronous  as  defined  above,  we  must 
restrict  to  be  of  the  form 


ri  =  n  Tq  (2.3) 

where  n  is  equally  likely  to  take  any  integer  value  on  the  closed  interval 
[0,  L-1].  We  designate  the  random  variable  for  the  cross-correlation  for 
integral  offset  by  c.  The  probability  mass  function  for  c  is  nonzero  only  for 
c  =  m  where  m  satisfies  -L  i  i  i  L  and  m  mod2  =  L  mod2.  For  an  arbitrary  timing 
offset,  can  be  represented  as 


t  =  (n  +  x)T  (2.4) 
i  c 

where  again  n  is  equally  likely  to  take  any  integer  value  on  the  closed 
interval  [0,  L-1]  and  x  is  uniformly  distributed  on  the  interval  (0,1).  In 
this  case,  the  code  cross-correlation  can  take  noninteger  values  due  to  partial 
agreements  and  disagreements  during  the  chip  intervals.  The  case  of  arbitrary 
timing  offsets  (to  include  nonintegral  offsets)  is  more  realistic  since,  in 
general,  it  is  neither  desirable  nor  feasible  to  impose  any  synchronization  on 
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the  system,  due  to  the  complexity  and  cost  which  would  be  involved.  Let  c^  be 
a  random  variable  corresponding  to  the  value  of  the  cross-correlation  when 
x  =  0,  and  c1  be  the  value  when  x  =  1.  For  given  values  of  cQ  and  c1 ,  the 
cross-correlation  varies  linearly  with  x  as  each  chip  comparison  transitions 
from  its  value  for  Cg  to  its  value  for  c^.  Consequently,  the  cross-correlation 
for  arbitrary  offset,  designated  c,  is  given  by 


c  =  (1  -  x)Cq  +  xc.  .  (2.5) 

When  Cg  =  c1 ,  c  takes  the  value  cQ.for  all  values  of  x.  Hence,  the  probability 
of  c  taking  the  value  k  is  equal  to  the  probability  that  both  cQ  and  e1  take 
the  value  k,  that  is, 


Pr{c  =  k}  =  Pr{Cg  =  k,  c1  =  k)  =  P£  k  (2.6) 

where  P'  is  the  joint  probability  mass  function  for  the  ordered  pair  (cn  =  r, 
c1  =  s).  To  obtain  the  density  for  c,  we  first  obtain  the  conditional  density 
of  c  given  cQ  and  c^.  Since  x  is  uniform  on  (0,1),  we  have 


pa(xlco*ei)  =  i 


|c1  -  c0|-1  ,  min(cQ,c1)  <  X  <  maxjcg.c^ 


(2.7) 


0  , 


otherwise  . 


The  density  is  then  obtained  by  taking  the  expectation  of  the  conditional 
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density.  Since  the  conditional  density  is  constant,  the  density  of  c  is 
piecewise  constant  with  discontinuities  occurring  whenever  X  =  k  for  -L  1  k  i  L 
and  k  mod2  =  L  mod2.  For  X  on  the  open  interval  (k,k+2)  where  k  is  even  for  L 
even  and  odd  for  L  odd,  we  have  (r  mod2  =  s  mod2  ■  L  mod2) 


k 

p-(X)  = 

r=-L 


L 


l 

s=k+2 


(a-r)"1 


r  ,s 


L 


l 

r=k+2 


k  , 
I  (r-3)"1 
s=-L 


r  ,s 


(2.8) 


where  any  open  sums  are  to  be  interpreted  as  zero.  This  equation  can  then  be 
manipulated  into  the  form  [Hanlon,  Peterson,  and  Gardner,  1980] 


P6(X) 


(L+k)/2  L 

l  .  I  . 

r-0  s=(L+k+2)/2 


(s-r)”1  P 


2r-L,2s-L 


where  (k  <  X  <  k+2,  -Llki  L-2,  k  mod2  =  L  mod2)  and 


(2.9) 


P 


r 


s,r 


(2.10) 


It  should  be  noted  that  for  the  Joint  probability  mass  function  P'  ,  the  order 

r ,  s 

of  r  and  s  is  important  since  r  designates  the  value  of  the  cross-correlation 

for  x  =  0  (cQ)  and  s  designates  the  value  for  x  =  1  (c^,  whereas,  in  Eq.(2.9), 

the  result  is  independent  of  the  order  since  P  is  the  average  of  P'  and 

r )S  r ,3 

P'  .  Consequently,  P„  can  be  considered  to  be  the  joint  probability  mass 
3|  r  r )  3 

function  of  c^  and  c^  without  regard  to  order. 


The  relations  shown  in  Eq.(2.6)  and  Eq.(2.9)  demonstrate  that  the 
distribution  of  c  is  dependent  only  on  the  joint  probability  mass  function  for 
adjacent  cross-correlation  values  under  the  condition  of  integral  offset.  This 
significantly  reduces  the  complexity  of  calculating  the  density  for  arbitrary 
offset  between  the  codes.  For  example,  consider  the  case  of  random  binary 
symmetric  codes,  that  is,  codes  for  which  the  value  of  each  chip  is  equally 
likely  to  be  either  +1  or  -1  and  for  which  each  chip  is  statistically 
independent  of  all  other  chips.  In  this  case,  for  r  mod2  »  s  mod2  »  L  mod2, 


P 

r,s 


2 


2L 

f 

L 

L 

L+r 

L+s 

2  J 

.  2  . 

and  from  Eq.(2.6)  and  Eq.(2.9),  we  have  for  the  mass  function 


(2.11) 


Pr{c  =  m}  =  2 


-2L 


L 

L+m 

2 


(2.12) 


where  (-Lim<  L,  m  mod2  =  L  mod2),  and  for  the  density 


Pg(X)  =  2 


■2L  (L41f/2 


L 

l 


r=0  s=(L+k+2)/2 


L 

f  > 

L 

r 

.3. 

(s-r)‘ 


(2.13) 


where  (k  <  X  <  k+2,  -L  i  k  i  L-2,  k  mod2  =  L  mod2).  Note  that  in  the  case  of 
random  codes,  the  magnitude  of  the  cross-correlation  can  be  as  large  as  L,  the 
code  length  (number  of  chips  per  data  bit). 
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For  a  general  set  of  periodic  codes,  determination  of  the  probability 
density  is  more  complicated.  Starting  with  Eq.(2.1)  and  noting  that  the  data 
bit  on  the  interfering  signal  has  a  constant  value,  say  b^,  on  the  interval 
(O.t^)  and  transitions  to  a  constant  value,  say  b^,  for  the  remainder  of  the 
integration  interval,  we  have 


c  ■-  61  V‘V  *  61  vv 


(2.1*) 


where  and  R^  are  partial  cross-correlation  functions  defined  by  [Pursley, 
1977] 


VV  *  T='  [ dt  *V-V  WV 


(2.15) 


Rji(Ti)  = 


„-1 


/  dt  aj(t-$j)  ai(t-$i-Ti) 


(2.16) 


Since  (b^,  b^)  can  assume  any  of  the  four  values  (♦1,+1)  with  equal  likelihood, 
then  for  a  given  offset  and  phase,  c  takes  on  the  four  values 


c  =  i  Rji(Ti)] 


(2.17) 


with  equal  likelihood.  It  is  readily  shown  by  a  change  of  variables  (see  also 
Massey  and  Uhran  [19753)  that  these  four  values  are  the  same  as 


26 


c  =  +[FiJ(Td-ri)  +  Rij(Td-Ti):i  •  (2.18) 

Then,  for  a  given  set  of  phases  for  and  a^,  varying  the  timing  offset  over 
all  possible  integral  values  generates  the  same  values  for  the 

cross-correlation,  probability  mass  function,  and  probability  density  whether 
a^  is  taken  to  be  the  intended  signal's  code  sequence  and  is  taken  to  be  the 
interfering  signal's  code  sequence  or  vice  versa. 

Note,  since  the  integration  interval  implied  in  Eq.(2.l4)  is  broken,  that 

the  set  of  values  for  the  partial  cross-correlation  functions  under  varying 

integral  offsets  changes  if  either  code  sequence  is  circularly  shifted  by  q 

chips  where  q  lies  on  the  interval  [0,L-1].  This  circular  shift  i3  equivalent 

to  the  code  sequence  being  started  with  a  different  phase  relative  to  its  data 

bit  sequence.  In  principle,  we  could  select  the  phase  for  each  code  to 

minimize  the  average  cross-correlation  moments  or  average  bit  error  rate. 

However,  for  long  codes  this  is  impractical,  and  we  do  not  assume  any 

particular  set  of  phases  for  the  codes,  that  is,  the  relative  phase  is  assumed 

to  take  on  any  of  the  L  possible  values  with  equal  likelihood.  We  can 

determine  the  values  of  P'  by  tabulating  the  relative  frequency  of  the 

r ,  s 

ordered  pair  (Cq.c^),  where  cQ  corresponds  to  t  3  nTc  and  corresponds  to 

=  (n+DT^,  under  exhaustive  enumeration  of  the  partial  cross-correlation 

functions  over  all  values  of  n  on  the  interval  [O.L-1];  over  all  code  pairs 

(a^.a^),  for  i  greater  than  j;  and  over  all  phases  of  the  codes.  The  symmetry 

of  Eq.(2.17)  and  Eq.(2.l8)  obviates  consideration  of  the  cases  where  i  is  less 

than  J.  P  is  then  the  average  of  P'  and  P'  .  The  values  of  P  are  then 
r ,s  r,3  s,r  r,s 

used  in  Eq.(2.6)  and  Eq.(2.9)  to  obtain  the  probability  density. 
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For  the  case  of  Gold  codes  of  length  L,  there  are  (L+2)  distinct  code 

sequences  which  can  be  generated  by  a  linear  shift  register  with  2m  stages, 

where  L  =  (2m  -  1).  For  example,  there  are  33  possible  codes  of  length  31 

which  can  be  generated  by  a  10  stage  shift  register.  Let  a^  represent  the 

(k) 

nominal  zero  phase  of  code  a^  relative  to  its  data  bit  and  a^  represent  a 

phase  shift  of  k  chips  relative  to  a^.  From  the  closure  property  of  linear 

shift  registers  (the  modulo  2  sum  of  two  sequences  from  a  linear  shift  register 

is  itself  a  sequence  from  that  shift  register)[Gold,  1964]  and  the  isomorphism 

between  the  additive  group  {0,1}  and  the  multiplicative  group  {-1,1},  it  is 

readily  shown  that  the  chip-by-chip  product  a^a^,  for  a^  not  equal  to 

a^m\  can  always  be  resolved  as  the  product  of  a^  and  a^3\  that  is,  the 

product  of  the  nominal  zero  phase  of  one  of  the  Gold  codes  and  some  other 

member  of  the  set  of  Gold  codes.  The  values  for  P  generated  by  the  code 

r  ,s 

combination  a^a^  then  replicate  the  set  of  values  generated  by  the  code 
combination  ap0^a£3^.  Consequently,  we  only  consider  phase  shifts  of  the 
interfering  signal's  code  sequence  rather  than  phase  shifting  both  the  desired 


and  interfering  signals'  code  sequences. 


Calculation  of  P  for  Gold  codes  then  involves 
r,s 


code  pairs,  L 


different  values  for  the  relative  timing  offset,  L  different  values  for  the 


relative  phase  of  a^,  and  4  different  data  bit  combinations.  Hence  there  are 
[2L  (L+1ML+2)]  values  for  the  cross-correlation.  Due  to  symmetry,  only  half  of 
these  values  need  be  calculated.  Consequently,  the  number  of  cross-correlation 
values  calculated  is  proportional  to  the  fourth  power  of  the  code  length. 
Figure  3  illustrates  the  results  for  a  set  of  Gold  codes  of  length  31  (m=5). 
Without  the  presence  of  the  data  bits,  the  cross-correlation  between  Gold  codes 
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of  length  L  is  three-valued:  (-2L (m+2)/2J  _  it  _i,  2L (m+2)/2 J  _  ^  [Gold> 
1964],  Figure  3a  shows  these  three  values  (-9,  -1,  7)  and  the  average  relative 
frequency  with  which  they  occur  for  integral  offsets  for  the  Gold  codes  of 
length  31.  The  results  shown  in  Fig.  3b  illustrate  the  effects  of  data  on  the 
cross-correlation  values  for  integral  offset.  Note  that  the  probability  mass 
function  is  symmetric  and  the  maximum  magnitude  of  the  cross-correlation  (21) 
extends  considerably  beyond  the  corresponding  value  without  data  present  (9). 
Figure  3c  shows  the  average  cross-correlation  probability  density  and 
probability  mass  function  for  arbitrary  offset.  As  mentioned  in  the 
development  for  arbitrary  code  sets,  the  density  is  piecewise  constant. 


2.3.  Moments 

Given  the  probability  density,  it  is  a  simple  matter  to  calculate  the 
moments  of  the  cross-correlation.  However,  if  we  desire  the  average 
performance,  the  number  of  code  correlations  needed  to  compute  the  density 
rapidly  becomes  prohibitive  as  the  code  length  increases.  Hence,  we  are 
motivated  to  develop  an  approximation  to  the  average  error  rate  which  can  be 
evaluated  without  use  of  the  probability  density.  One  approach  is  to  expand  the 
conditional  probability  of  error  given  c  in  a  power  series  of  the  form 


Pr{E|c}  *  l  A  (c/L)2k  (2.19) 

k  * 

where  the  A^  are  the  coefficients  of  the  expansion.  Provided  such  a  series 
exists  and  that  it  converges  sufficiently  rapidly,  the  expectation  of  the  first 


Figure  3.  Probability  mass  functions  and  density  for  cross-correlation  of 
Gold  codes  of  length  31:  (a)  relative  frequency  of  values  for 

integral  offset,  without  data:  (b)  mass  function  for  integral 
offset,  with  data;  (c)  mass  function  and  density  for  arbitrary 
offset,  with  data. 


few  terms  provides  an  accurate  estimate  of  the  error  rate. 


This  would  limit 


the  problem  to  determining  the  first  few  moments  of  the  cross-correlation  or, 
more  practically,  an  estimate  of  these  moments.  This  approach  has  been  shown 
to  provide  accurate  estimates  of  the  probability  of  error  for  direct  sequence 
systems  operating  over  nonselective  slow  fading  channels  [Gardner  and  Orr, 
1979],  [Hanlon  and  Gardner,  1979  and  1980]. 


It  has  been  previously  conjectured  that  the  low-order  moments  of  the 
average  cross-correlation  in  the  presence  of  data  are  close  to  those  for  random 
code  sequences  of  the  same  length.  It  is  shown  in  Appendix  D  that  the 
cross-correlation  moments  of  order  2n  for  random  code  sequences  can  be 
expressed  as  a  polynomial  of  degree  n  in  the  code  length.  Table  D.2  lists  the 
coefficients  of  these  polynomials  for  the  first  five  even  moments.  For  large 
values  of  L  (long  codes),  the  moments  are  shown  to  be  approximately  given  by 


,~2n. 
<c  > 


(2n  -  1 ) ! !  L,n  y  fn)  (2n]  1 

2n+l  mtQ  (mj  (2;nJ 


[m2+(n-m)2-n] /3L}  . 


(2.20) 


This  expression  is  exact  for  the  second  and  the  fourth  moments  which  are, 

2 

respectively,  (2L/3)  and  (7L  /5  -  9L/5).  These  values  differ  significantly 

2 

from  the  corresponding  moments  for  integral  offset,  which  are  L  and  (3L  -  2L) 

(Table  D.1).  Consequently,  when  analyzing  the  performance  of  nonsynchronous 
systems,  estimates  of  the  probability  of  error  or  related  functions  (e.g., 
signal-to-noise  ratio)  which  make  use  of  the  cross-correlation  moments,  must 


consider  arbitrary  offsets  rather  than  just  integral  offsets. 
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Tables  1  and  2  show  the  low-order  moments  of  random  codes  and  Gold  codes 
for  code  lengths  of  31  and  127,  respectively.  As  anticipated,  the 
corresponding  moments  are  fairly  close,  particularly  the  second  moments. 
However,  as  the  order  of  the  moments  increases,  the  values  begin  to  differ 
significantly.  The  randomizing  effects  of  the  interfering  signal's  data 
sequence  and  the  relative  timing  offset  cause  the  density  of  the 
cross-correlation  for  Gold  codes  in  the  presence  of  data  to  take  a  similar  form 
to  that  of  random  code  sequences  of  the  same  length.  One  significant 
difference  is  the  maximum  value  of  the  cross-correlation  magnitude.  The  effect 
of  this  difference  on  the  second  moment  for  Gold  codes  relative  to  random  codes 
is  small  due  to  the  low  probability  that  the  cross-correlation  for  random  codes 
exceeds  the  maximum  value  for  Gold  codes.  However,  as  the  order  of  the  moments 
increases,  the  outlying  values  are  weighted  more  in  the  calculation  of  the 
moments  and  the  difference  in  values  for  the  moments  is  progressively  more 
pronounced. 

To  facilitate  development  and  use  of  the  empirical  estimates,  it  is 
desirable  to  consider  a  continuous  density.  If  we  develop  a  close 
approximation  to  the  low-order  Gold  code  moments  by  use  of  the  moments  for  a 
continuous  density,  we  would,  in  principle,  be  able  to  extrapolate  the 
approximation  to  provide  estimates  of  the  low-order  absolute  moments  or 
low-order  nonintegral  moments  of  the  Gold  codes.  These  generalized  moments  are 
of  potential  use  when  the  series  expansion  of  the  conditional  probability  of 
error  given  c  has  a  more  complicated  dependence  on  c  than  that  shown  in 
Eq . ( 2 . 1 9 ) . 


TABLE  1.  CROSS-CORRELATION  MOMENTS  <c2n>  (L  =  30 

n  Gold  Codes  Random  Codes  Gamma  Density 

1  2.05  E  1  2.07  E  1  2.07  E  1 

2  1.12  E  3  1.32  E  3  1.15  E  3 

3  9.24  E  4  1.45  E  5  9.95  E  4 

4  1.08  E  7  2.27  E  7  1.15  E  7 

5  1.73  E  9  4.64  E  9  1.*5  E  9 


TABLE  2.  CROSS-CORRELATION  MOMENTS  <c2n>  (L  =  127) 

n  Gold  Codes  Random  Codes  Gamma  Density 

1  8.75  E  1  8.47  E  1  8.47  E  1 

2  1.97  E  4  2.25  E  4  1.93  E  4 

3  6.58  E  6  1.04  E  7  6.84  E  6 

4  2.96  E  9  6.99  E  9  3.24  E  9 

5  2.04  E12  6.26  E12  1.90  E12 
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The  symmetry  and  general  shape  of  the  density  for  c  (Fig.  3)  suggest  a 
Gaussian  approximation.  The  main  disadvantage  of  using  a  Gaussian  density  is 
that  the  range  of  the  density  is  unbounded  (tails  extend  indefinitely).  This 
tends  to  exacerbate  the  effect  of  the  difference  in  the  magnitude  of  the 
maximum  cross-correlation  values.  However,  the  density  falls  off  fairly 
rapidly  and  would  provide  reasonable  approximations  to  the  low-order  moments, 
which  are  of  prime  concern  to  our  application.  The  effect  of  the  tails  can  be 
reduced  if  the  Gaussian  density  is  generalized  slightly  to  increase  the  rate  at 
which  the  tails  decay.  One  generalization  of  the  Gaussian  density  is  the 
two-sided  generalized  gamma  density.  The  form  of  the  generalized  gamma  density 
is  given  by 

P  ( x )  =  A  |x  I01-"1  exp(-e|x[v)  (2.21) 

where  for  the  Gaussian  density,  a  equals  one  and  v  equals  two.  By  keeping  a 
equal  to  one  and  setting  the  second  moment  equal  to  that  of  the  random  codes, 
i.e.,  2L/3,  the  moments  for  the  gamma  density  are  then  given  by 

<x2n>  =  [r(^2±l)/r(i)]  [2r(l)/3r(l)]n  Ln  .  (2.22) 

v  v  v  v 

To  develop  the  empirical  estimates,  the  value  of  v  is  then  chosen  to 
approximately  match  the  first  five  even  moments  of  the  gamma  density  to  the 
corresponding  moments  of  the  cross-correlation  for  Gold  codes  of  length  31. 
The  appropriate  value  of  v  is  about  2.4.  It  is  expected  that  this  choice  of 
v  would  also  result  in  a  good  estimate  of  the  low-order  moments  for  Goli  ^odes 
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of  longer  lengths  due  to  (1)  anticipated  close  agreement  of  the  second  moments, 
regardless  of  code  length,  and  (2)  the  dependence  on  code  length  for  these 
estimates  is  the  same  as  in  the  single-term  approximation  for  the  cross¬ 
correlation  moments  for  random  codes.  Table  1  shows  the  moments  for  the 
average  cross-correlation  for  Gold  codes  of  length  31 ,  the  corresponding 
moments  for  random  sequences  of  the  same  length  and  moments  of  the  generalized 
gamma  density  (v  =  2.4).  The  moments  of  the  generalized  gamma  density  are  seen 
to  agree  very  well  with  the  moments  of  the  Gold  codes. 

For  Gold  codes  of  length  127,  the  moments  were  calculated  using  a  subset  of 
5  out  of  the  possible  129  Gold  codes  (permitting  10  pairwise  combinations  out 
of  a  possible  8256  combinations) .  We  are  still  concerned  primarily  with  the 
average  over  the  entire  code  set;  however,  to  reduce  the  calculations,  we  have 
used  this  subset  to  approximate  the  desired  results.  Table  2  shows  the 
corresponding  moments  for  this  subset.  The  same  general  pattern  appears  as  in 
the  case  of  Gold  codes  of  length  31.  The  agreement  in  the  second  moments  is 
very  good  considering  the  small  sample  used. 

The  random  code  moments  are  of  particular  interest  in  related  applications. 
In  order  to  take  full  advantage  of  the  potential  benefits  of  spread  spectrum 
systems,  very  long  code  sequences  are  desirable.  These  long  code  sequences 
enhance  the  privacy  and  anti jam  characteristics.  For  a  given  type  of  code,  the 
longer  the  code  sequence,  the  harder  it  is  to  effectively  monitor  or  jam  the 
communications  system  on  which  it  is  used.  However,  there  are  practical  limits 
(technical  and  economic)  to  arbitrarily  expanding  the  transmitted  bandwidth 
[Dixon,  1976a],  In  some  applications  then,  the  period  of  the  spreading  codes 


is  chosen  to  be  much  greater  than  the  data  bit  interval.  This  provides  for 


both  the  maximum  feasible  degree  of  spectrum  spreading  and  very  long  code 
sequences.  Let  L  be  the  period  of  the  spreading  codes  and  L'  be  the  number  of 
chips  per  data  bit,  where  L  is  much  greater  than  L'.  Then  if  the  spreading 
codes  exhibit  a  high  degree  of  randomness  [Golomb,  1967]  over  a  data  bit 
interval,  such  as  in  the  case  of  cryptographic  codes,  we  would  expect  that  the 
cross-correlation  moments  would  be  very  closely  approximated  by  the  moments  of 
random  codes  of  length  L'. 
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Chapter  3 

PROBABILITY  OF  ERROR 

3.1.  Single  Interferer  Without  Additive  Noise 

For  the  case  of  a  single  interferer  where  the  additive  noise  is  negligible, 
the  probability  of  error  represents  an  irreducible  error  rate  which  bounds  the 
multiple-access  system  performance.  Again,  the  restriction  of  consideration  to 
a  single  interferer  does  not  imply  that  there  is  only  one  other  system  user, 
but  rather  that  only  one  of  the  other  system  users  is  simultaneously 
transmitting.  In  this  case,  from  Eq.(1.12),  the  probability  of  error  is  given 
by 


?r{E}  =  £  Pr{af2  >  f 1 } 


(3.1) 


where 


a  =  (Pg/P,) 


C  COS  0 
n  1 


5  =  (c/L)  ,  0=0,,. 

n  2 


(3.2) 


It  should  be  noted  that  f 1 ,  fg,  cn,  and  0  are  all  random.  To  facilitate  the 
evaluation  of  this  expression,  the  conditional  probability  of  error  given  the 
factor  a  is  first  determined. 
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It  is  readily  shown  that  if  f  has  a  generalized  Rayleigh  distribution,  that 

is,  f  ^  R(m,g,ft),  then  the  product  af  also  has  a  generalized  Rayleigh 

2 

distribution:  af  ^  R(m,aS,a  fi).  From  Appendix  B,  Eq.(B.26),  we  have 


00 

Pr{E  j a }  =  ^  (1-p2)m  expt-mg2/^)  £  '(a2^/^)1 


(-1)k  Ok(m,D2/B1,D2/B2,p2) 


(3.3) 


where  Gk  is  given  by  Eq.(B.27)  and  B1 ,  B2,  D1 ,  and  D2  are  given  by  Eq.(B.21). 
It  should  be  noted  that  the  function  Gk  has  no  dependence  on  the  factor  a.  The 
probability  of  error  is  then  obtained  by  taking  the  expectation  of  the 
conditional  probability  of  error  with  respect  to  both  &n  and  0  which  are 
independent  of  each  other.  From  the  definition  of  the  beta  function  and  the 
assumption  that  0  is  uniform  on  the  interval  (0,2-it),  we  readily  obtain 


<  [ cos  0|2k+2m>  =  B(1/2,k+m+1/2)  /  it  =  m”1/2  --(m+1-)  k 


By  applying  this  result  to  Eq.(3.3),  we  have  [Hanlon  and  Gardner,  1980] 
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The  bivariate  density  given  by  Eq.(B.l8)  implicitly  assumes  that  the  m 
factors  are  the  same  for  both  fading  intensities.  If  we  assume  that  the  fading 
intensities  are  independent,  so  that  the  joint  density  is  given  by  the  product 
of  the  marginal  densities,  we  can  readily  extend  the  results  to  the  case  where 
the  m  factors  are  different  (Eq.(B.II))  [Hanlon  and  Gardner,  1980] 
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(3.6) 


The  practical  usefulness  of  Eq.(3-5)  or  Eq . C 3 . 6 )  depends  on  the  series 
converging  sufficiently  rapidly  to  enable  close  approximation  of  the  average 
probability  of  error  by  one  or  two  terms  of  the  series.  Since  the 
cross-correlation  moments  for  the  spreading  codes  are  typically  very  small  for 
well-designed  DS/SSMA  systems,  rapid  convergence  is  quite  feasible  provided  the 
ratio  (m1P2f22/m2P1S11 )  is  not  large.  Table  3  lists  the  single-term 
approximations  for  the  more  important  special  cases  of  the  generalized  Rayleigh 
distribution. 

Numerical  evaluation  of  the  average  probability  of  error  requires 
determination  of  the  absolute  moments  of  the  normalized  code  cross-correlation. 


This  can  be  very  difficult,  particularly  since  for  generalized  Rayleigh  or 
Nakagami  fading,  we  require  nonintegral  absolute  moments.  In  general,  we  would 
need  the  probability  density  of  the  code  cross-correlation;  however,  for  any 
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TABLE  3 

SINGLE-TERM  APPROXIMATIONS  FOR  THE  PROBABILITY  OF  ERROR  FOR  SPECIAL  CASES  OF 
GENERALIZED  RAYLEIGH  FADING  FOR  A  SINGLE  INTERFERER 

f  2 

t(1-P2)/U]  exp( -8^/12  )  (p  a  /p  a  )  <C2>  1  - 2 

11  2  2  11  n  nffi-d-a2) 

l  1  1 

where  (^  =  8 2  “  p  ^  fi]^2  n2^2 

NAKAGAMI:  f2”;1)2  [ ( 1-P2 )/M ]m  (Pg^/P^)®  <!cn|2m> 

(m1sm2=m,  81=82=0) 


RICIAN: 
(nj1=m2=1 ) 


RAYLEIGH:  C(1-P2)/4]  (P,n_/P  ,0,)  <c2> 

2  2  1  1  n 

(m1=m2  =  1,  8^82=0) 


INDEPENDENT  FADING: 
(p=0) 


2m, -1 

."l 


-m  -m  m-, 

exp(-m1  8^/^)  2  (m1P2^2/m2Pini^ 


(m2_1)  A  I 

L  )  <|8n 

m1  2  2  2  n 


Note:  for  Gold  codes,  <c2>  =  2/3L  and  <  |c  |2m>  is  approximately  given  by 
Eq.(2.22),  after  normalization. 


reasonably  long  code,  determination  of  the  probability  density  is 
computationally  impractical.  Therefore,  estimates  of  the  moments  are  desired. 
For  the  case  of  Gold  codes,  an  empirical  estimate  of  the  low-order  even  moments 
of  the  cross-correlation  is  given  by  Eq.(2.22).  By  straightforward 
extrapolation  of  this  result  to  nonintegral  absolute  moments,  we  can  readily 
evaluate  the  single-  or  two- term  estimates  of  the  average  probability  of  error, 
provided  the  parameter  m  is  not  too  large,  i.e.,  provided  the  order  of  the 
moments  required  is  not  too  large. 

Figures  4  through  11  are  plots  of  the  average  probability  of  error,  the 

single-term  estimates  using  the  approximate  moments  given  by  Eq.(2.22),  and  the 

corresponding  aproximation  error  for  several  examples  of  Rician  and  Nakagami 

2  2 

fading.  The  curves  are  drawn  for  the  case  of  equal  power  (P1 <f1 >=P2<f2>) , 
equal  fading  statistics  (m1=m2,  6^s0o,  ,  and  a  subset  of  Gold  codes  of 

length  127.  The  plots  of  the  exact  values  for  the  probability  of  error  were 
obtained  by  evaluating  the  moments  from  the  actual  probability  density  for  the 
code  cross-correlation  taken  over  a  subset  of  5  of  the  129  Gold  codes  of  length 
127  (permitting  10  pairwise  combinations  out  of  a  possible  8256  combinations). 
The  probability  density  is  determined  by  tabulating  the  relative  frequency  of 
cross-correlation  values  for  all  possible  integer  offsets  (T1=nTc,  0^.n<L)  for 
each  of  the  code  pairs,  for  all  possible  data  sequences  on  the  interfering 
signal  (b^),  and  for  all  initial  phases  of  the  interfering  signal's  spreading 
code  (circular  shifts  of  the  code  within  the  data  bit  interval)  (Chapter  2], 
[Hanlon,  Peterson,  and  Gardner,  1980].  This  was  practicable  due  to  the  short 
code  length  and  the  limited  code  set.  The  plots  show  that  the  single-term 
estimates  are  quite  good  over  a  wide  range  of  the  channel  fading  parameters. 


CORRELATION  COEFFICIENT, rho 


CORRELATION  COEFFICIENT,  rho 


RICIAN  SPECULAR  POWER 


MAGNITUDE  OF  CORRELATION  COEFFICIENT 


MAGNITUDE  OF  CORRELATION  COEFFICIENT 


NAKAGAMI  M- DISTRIBUTED  FADING  FACTOR 
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NAKAGAMI  M -DISTRIBUTED  FADING  FACTOR 
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The  Rician  results  are  given  for  specular  power  ratios  of  0,  50,  67,  and  75 

2 

percent,  which  correspond  to  power  factors  [R/( 1  — R )  =  $  /Q]  of  0,  1,  2,  and  3- 
As  the  correlation  of  the  underlying  Gaussians  increases  (  p  increases) ,  the 
correlation  between  f^  and  f^  increases  and  the  probability  of  error  drops 
towards  zero.  As  the  fading  factors  more  closely  agree  statistically,  the 
signal  and  the  interference  tend  to  maintain  their  relative  strengths  at  the 
input  to  the  receiver.  As  a  result,  it  is  less  likely  that  the  sum  of  the 
signal  and  the  interference  will  result  in  a  sign  reversal  in  the  detector  and 
cause  a  bit  error  [Hanlon  and  Gardner,  1979].  This  is  also  true  for  Nakagami 
fading  and,  in  general,  the  import  of  the  fading  is  significantly  reduced  for 
highly  correlated  channels.  As  the  Rician  specular  power  ratio  increases 
(fading  decreases),  the  probability  of  error  decreases;  however,  the 
approximation  error  tends  to  increase.  Use  of  more  than  one  term  of  the  series 
in  the  approximation  would  result  in  a  wider  range  of  the  parameters  over  which 
close  agreement  is  attained. 

The  Nakagami  results  are  very  similar  to  those  for  Rician  fading.  The 
error  rate  decreases  as  the  fading  decreases  (m  increases)  and  the 
approximation  error  tends  to  correspondingly  increase.  The  empirical  estimates 
of  the  code  cross-correlation  moments  can  be  used  to  calculate  additional  terms 
of  the  series  to  improve  the  accuracy,  provided  m  is  not  made  larger  than  about 
4.  For  larger  values  of  m,  the  order  of  the  moments  required  are  too  large  to 
obtain  good  estimates  of  the  moments  from  the  empirical  moments  [Hanlon, 
Peterson,  and  Gardner,  1980], 


all 
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The  average  probability  of  error  for  Rician  fading  is  dominated  by  the 

depth  of  fading  on  the  intended  signal  as  evidenced  by  the  exponential  factor, 
2 

exp(-81/01).  a  similar  dependence  for  Nakagami  fading  is  reflected  in  the 
exponential  dependence  on  m  (order  of  the  moments  of  the  code 

cross-correlation) .  Since  we  have  implicitly  assumed  that  the  m  factor  is  the 
same  for  both  fading  intensities,  the  distinction  is  difficult  to  see.  For  the 
case  of  independent  Nakagami  fading,  Eq . ( 3 - 6 )  clearly  shows  that  the 
exponential  dependence  is  on  (fading  on  intended  signal)  rather  than 

(fading  on  interfering  signal). 

While  the  performance  for  both  the  Rician  and  Nakagami  fading  cases 
demonstrates  high  dependence  on  the  likelihood  of  deep  fades  on  the  intended 
signal,  the  performance  shows  significantly  different  dependence  on  the  code 
length  in  the  two  cases.  This  is  readily  seen  by  comparing  the  single-term 

a2 

approximations.  For  Rician  fading,  the  approximation  is  proportional  to  <cn> 
and,  hence,  from  Eq.(2.22),  L_1 ,  whereas,  the  single-term  approximation  for  the 
Nakagami  case  is  proportional  to  <|cn|^m>,  or  equivalently,  L-m.  The  basis  for 
this  difference  can  be  seen  by  looking  at  the  case  of  independent  channels. 
For  independent  channels,  the  error  rate  can  be  alternatively  represented  as 
[Appendix  C],  [Hanlon  and  Gardner,  19791 


PrtEl  =  ~  l  <|c!k>  <!cos  ©|k>  (P./P.)k/2 


<fk>  pi.k_1)(0) 
d  1 1 


(3.7) 


lr 

where  <f^>  is  the  k-th  moment  of  the  fading  intensity  on  the  interfering  signal 

(k-i ) 

and  p  (0)  is  the  (k-1)  derivative  of  the  probability  density  for  the  fading 
*  ■! 
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intensity  of  the  intended  signal  evaluated  at  zero.  This  expression  requires 
the  derivatives  at  the  origin  to  be  finite,  which,  in  general,  is  not  the  case 
for  the  generalized  Rayleigh  distribution.  By  restricting  our  interest  to 
integral  values  of  (23^),  we  can  see  by  expanding  the  probability  density  for 
f1  in  a  power  series  that  the  first  nonzero  derivative  at  the  origin  is  for  k 
equal  to  2my  The  larger  the  value  of  (the  shallower  the  fading  on  the 
intended  signal),  the  higher  the  order  of  the  code  cross-correlation  moment 
which  first  contributes  to  the  series.  This  is  exactly  the  effect  seen  in 
Eq.(3.6),  which  can  also  be  derived  from  Eq.(3.7)  when  (23^)  is  an  integer. 
From  Eq.(2.22)  we  see  that  the  first  term  varies  as  L-m.  For  Rician  fading,  m 
is  equal  to  one  and  the  first  term  is  always  proportional  to  the  normalized 
second  moment,  which  is  to  say,  inversely  proportional  to  the  code  length. 
Consequently,  the  degree  of  improvement  in  performance  to  be  expected  from 
increasing  the  code  length  is  highly  dependent  on  the  form  of  the  distribution 
of  the  intended  signal  fading  intensity,  particularly  the  behaviour 
(derivatives)  of  the  distribution  near  the  origin  (region  of  deep  fades).  As  a 
result,  we  cannot,  in  general,  approximate  a  given  fading  distribution  with 
another  distribution.  For  example,  Nakagami  gives  the  transformation  for 
approximating  a  Rician  distribution  by  a  Nakagami  distribution  as  [Nakagami, 
I960],  [Nakagami,  et  al.,  1957] 


=  nR  ♦  b2  ,  m  =  [i  -  (s4/  n*)]"1  • 


(3.8) 


This  transformation  is  obtained  by  equating  the  second  and  fourth  moments 
(Eq.(B.4)).  In  terms  of  the  specular  power  ratio  (Eq.(B.6)),  this  implies 
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m  =  (1  -  R2)-1  .  (3.9) 

Under  this  transformation,  specular  power  ratios  of  0,  50,  67,  and  75  percent 
correspond  to  m  values  of,  respectively,  1,  1.33,  1.8,  and  2.29.  It  is  clear 
from  Fig.  4  and  Fig.  8  that  this  transformation  does  not  provide  comparable 
performance  measures  in  terms  of  bit  error  rates.  Likewise,  any  transformation 
which  depends  primarily  on  equating  moments  and,  therefore,  on  the  overall 
behaviour  of  the  distribution,  generally  provides  poor  agreement  since  the 
error  rate  is  determined  primarily  by  the  behaviour  of  the  distribution  in  the 
region  of  deep  fades  [Barrow,  1962].  Regardless  of  the  form  of  the 
distribution,  relatively  long  code  sequences  are  required  to  maintain  low  error 
rates  unless  the  fading  on  the  intended  signal  is  very  shallow  or  the  fading 
correlation  approaches  unity. 

3.2.  Multiple  Independent  Interferes  With  Additive  Noise 

For  the  case  of  several  simultaneous  interferers,  the  individual 
interference  terms  in  Eq.(1.9)  are  independent  of  each  other  due  to  the 
presence  of  their  respective  data  bits  and  timing  offsets.  The  probability 
density  for  each  term  will  be  symmetric  due  to  the  cross-correlation  factors 
and  phases.  These  factors  also  tend  to  scale  the  individual  terms  towards  zero 
(  |cn  |  «  1,  | cos  0|  £  1),  i.e.,  towards  their  means.  Provided  the  P  are  on  the 
same  order  (which  will  be  assumed),  then  none  of  the  individual  terms  dominates 
the  series.  Under  these  conditions,  the  composite  interference  at  the  output 
of  the  correlator  (Z)  will  be  asymptotically  normal  with  zero  mean  and  variance 
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Var(Z)  =  Var(N)  +  Var(I)  =  (n„T/l|) 
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(3.10) 


The  probability  of  error  is  (Eq .(1.12)) 


PrtE}  =  j  Pr{  |Z  |  >  S}  . 


(3.11) 


Since  Z  i3  assumed  to  be  zero-mean  Gaussian,  |Z|  is  one-sided  Gaussian,  that 
is,  generalized  Rayleigh:  [Z|  'v  R[1/2,0,Var(Z)] .  Evaluation  of  the  probability 
of  error  then  reduces  to  determination  of  the  probability  that  one  generalized 
Rayleigh  variate  is  greater  than  another  when  they  are  independent.  From 
Eq. (B.16 ) ,  we  have 
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where,  for  M  system  users, 
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and  SNRq  is  the  signal-to-Gaussian-noise  ratio  defined  by 


(3.14) 


3NRq  =  [P^2  /  2  Var(N)]  =  2P1Td/  nQ  . 


(3.15) 
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The  usefulness  of  the  series  in  Eq.(3.12)  is  restricted  to  values  of  y  which 
are  less  than  one.  For  values  of  y  greater  than  one,  the  series  does  not 
converge.  There  are  two  approaches  to  circumvent  this  problem.  First,  we  can 
evaluate  the  complement  of  the  probability  to  obtain  an  expression  for  the 
probability  of  error  when  y  is  greater  than  one,  that  is, 


Pr{E}  =  [1  -  Pr{S  >  |Z [}]  /  2 


=  [1 


-1/2 
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Alternatively,  we  can  use  the  change  of  variables 


t  =  y/( 1  +  y)  . 


(3.17) 


By  using  the  binomial  theorem  and  Eq.(A.27),  we  then  obtain 
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Since  t  is  always  less  than  one  and  less  than  y,  the  convergence  properties  of 
this  expansion  are  much  better.  For  small  values  of  t,  i.e.,  t  much  less  than 
one,  the  average  probability  of  error  is  very  closely  approximated  by 


If  the  signal-to-Gausaian-noise  ratio  is  made  arbitrarily  large,  we  have 

lim  y  =  cun.  (M-1)  <c2>  (1  +  S2/n  )  .  (3.20) 

SNR  n 

G 

^2 

For  reasonably  long  code  sequences,  <cn>  is  very  small  and  y  tends  to  be  small. 
In  this  case,  t  and  y  are  approximately  equal  and  Eq.(3.19)  very  closely 
approximates  the  average  probability  of  error.  Since  we  are  considering  an 
arbitrarily  large  signal-to-Gaussian-noise  ratio,  this  represents  an 
irreducible  error  rate  due  to  cochannel  interference  and  fading.  Table  4  lists 
the  estimates  of  the  irreducible  error  rates  for  the  various  special  cases. 

We  have  analyzed  the  multiple  interferer  case  by  assuming  that  the 
composite  interference  is  asymptotically  normal  (Gaussian)  with  zero  mear.  As 
a  result,  the  interference  was  specified  by  it3  variance  (power).  We  would  not 
typically  expect  the  Gaussian  assumption  to  hold  if  the  number  of  interferers 
was  very  small  (say  less  than  five).  This  can  be  tested  by  looking  at  the 
single-term  estimates  for  the  irreducible  error  rate,  evaluating  them  for  a 
single  interferer  (M=2),  and  comparing  the  estimates  with  the  results  of 
Section  3.1  for  a  single  interferer.  The  single-term  estimate  for  the 
irreducible  error  rate  for  multiple  users  when  M  equals  two  is  given  by 
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TABLE  4 

SINGLE-TERM  APPROXIMATIONS  FOR  THE  IRREDUCIBLE  ERROR  RATE  FOR  MULTIPLE 
INDEPENDENT  USERS 
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Comparison  with  the  single-term  estimate  for  independent  fading  given  in  Table 

3  indicates  that  there  can  be  a  significant  difference.  However,  for  Rician 

fading  on  both  channels  (m^m^l),  the  result  given  by  Table  3  is  identical  to 

that  given  by  Eq.(3.21)  for  Rician  fading  on  the  intended  signal  (m. =1).  This 

comparison  extends  to  Rayleigh  fading  when  8^  ia  set  equal  to  zero  in  the 

Rician  estimates.  This  agreement  arises  because  the  single-term  single- 

interferer  estimates  in  both  cases  depend  only  on  the  interference  power.  In 

these  cases,  the  Gaussian  assumption  for  the  interference  i3  valid  even  when 

the  number  of  interferers  is  minimal.  However,  comparison  of  the  Nakagami 

results  (s.|  =  B2=0)  does  not  provide  the  same  agreement  unless  we  look  at  the 

special  case  of  Rayleigh  fading.  The  difference  arises  because  the  Nakagami 

multiple  user  result  is  dependent  on  the  m-th  power  of  the  interference  power, 

2  tn 

that  is,  proportional  to  (<cn>  );  whereas,  the  single  interferer  result 
(Section  3.1)  depends  on  the  (2m)-th  absolute  moment. 


Although  the  expression  for  the  probability  of  error  does  not  appreciably 
simplify  for  the  Rician  case  (m^=1),  for  Nakagami  fading  statistics  for  the 
intended  signal  (8^=0),  Eq.(3.12),  Eq.C3.l6),  and  Eq .(3.18)  reduce  to 


Pr{E} 


m^-1/2 


2F1(1/2,m1  ;m1  +  1;t)  =  j  I^m^l/2) 


(3-22) 


where  ^  is  a  hypergeometric  function  and  Ifc  is  a  normalized  incomplete  beta 
function.  This  expression  is  equivalent  to  a  result  given  by  Barrow  [1962]. 
It  is  readily  shown  that 


.1 
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2F1 (1/2,m;m+1;t)  =  (1-t) 


-1/2 


2(m+1 ) 


2F1 (3/2, m+1 ; m+2 ; t ) 


(3.23) 


By  truncating  both  hypergeometric  aeries  in  this  relation,  we  obtain  the  bounds 


[  1 


mt 


2(m+1) 


]  < 


^  ( 1  /2 ,  m ;  m+ 1 ;  t ) 


<.  [  (1-t) 


-1/2 


2  ( m+ 1 ) 


(3.24) 


Extension  of  these  bounds  to  the  probability  of  error  for  Nakagami  fading  i3 
straightforward  and  provides  the  estimate  (for  small  t) : 


Pr{E} 


[Vi/2 

l  "1 


[  1  +  (1-t) 


-1/2 


t(m1-1 ) 
2(m1+l)  ] 


(3.25) 


For  the  case  of  Rayleigh  statistics  (m^l,  8^0),  Eq.(3.22)  reduces  further  to 
the  well-known  result: 


Pr{E)  =  2  [1  "  =  g  [1  -  (1+y)"1/2l  •  (3.26) 


Figures  12  and  13  are  plots  of  the  average  probability  of  error  versus  the 
signal-to-Gaussian-noise  ratio  for  the  multiple  user  case  for  Rician  and 
Nakagami  (m-distributed)  fading,  respectively.  Both  sets  of  curves  are  drawn 
for  equal  effective  power  (P  =2^,  mean  square  values  normalized  to  one 
(<f?>  =  <f^>  =1),  and  <c^>  =  2/3L.  The  figures  clearly  demonstrate  that  the 
average  probability  of  error  saturates  at  an  irreducible  error  rate  for 


59 


O 

<r 


O 


o 

OJ 


O 


o 


o 


o  —  cj  m  «r  » 

o  g  o  'g  c 


u  n 

0  v 

^  p 


a  —  • 
^  cc 
o  u 

0)  W  OJ 

CO  u 

-4  <U 
C  •  iu 


CJ 

5} 


l 

1 


aoaaa  jo  Ainiavaoad 


SIGN AL-TO-GAUSSI AN- NOISE  RATIO  (dB) 


Probability  of  error  versus  s  I  f;na  l-to-flaiiss  lan-noi  se  ratio  for 
Nak;if',ami  fading  with  varying  values  of  llie  m-cll  str  i  bated  fad  in] 
factor.  Solid  lines  for  10  Interferers  anil  dashed  lines  for 
20  interferers. 


arbitrarily  large  signal-to-Gaussian-noise  ratios.  This  effect  is  due  to  the 
combination  of  cochannel  interference  and  fading  on  the  intended  signal.  The 
deeper  the  fading  on  the  intended  signal  or  the  greater  the  number  of 
interferers,  the  more  likely  that  the  intended  signal  component  falls  below  the 
level  of  the  cochannel  interference  component  in  the  receiver,  thereby 
resulting  in  an  error. 


Chapter  4 

SUMMARY,  CONCLUSIONS,  AND  RECOMENDATIONS 


4.1.  Summary 

This  study  has  investigated  the  error  performance  of  direct  sequence  spread 
spectrum  multiple-access  communications  systems  on  nonselective  fading 
channels.  The  receiver  model  used  throughout  the  analysis  is  a  coherent 
correlation  receiver.  The  receiver  correlates  the  received  signal  against  a 
synchronized  copy  of  the  intended  signal's  spreading  code  sequence.  The 
channel  model  assumes  slow,  nonselective  fading  where  the  fading  intensities 
are  distributed  as  generalized  Rayleigh  variates.  Results  are  obtained  for  the 
special  cases  of  a  single  interferer  without  additive  noise  and  multiple 
independent  interferers  with  additive  noise.  The  results  are  developed  as  a 
series  expansion  in  the  code  cross-correlation  moments  (single  interferer)  or 
as  a  series  expansion  in  the  powers  of  the  code  cross-correlation  second  moment 
(multiple  interferers).  To  develop  simple  estimates  of  the  probability  of 
error,  empirical  estimates  of  the  low-order  cross-correlation  moments  of  the 
spreading  codes  (taken  to  be  Gold  codes)  are  developed  from  a  gamma  density 
approximation  based  on  the  cross-correlation  moments  for  both  Gold  codes  and 
random  codes.  Numerical  results  for  the  probability  of  error  are  obtained  with 
these  estimates  and  compared  with  the  results  based  on  exact  moments  developed 
from  the  code  cross-correlation  probability  density. 


The  analysis  developed  in  this  paper  has  shown  that  the  method  of  series 
expansion  in  the  code  cross-correlation  moments  provides  a  tractable  solution 
to  evaluation  of  the  error  performance  of  direct  sequence  multiple-access 
communications  systems  operating  over  nonselective  fading  channels.  In 
particular,  truncation  of  the  series  provides  simple  single- term  estimates 
which  are  accurate  over  a  wide  range  of  the  channel  fading  parameters.  These 
results  show  that  the  error  rate  is  highly  dependent  on  the  form  of  the 
intended  signal's  fading  distribution  in  the  region  of  deep  fades. 
Consequently,  approximation  of  a  given  distribution  by  another  distribution 
(e.g.,  approximation  of  a  Rician  distribution  with  a  Nakagami  distribution)  is 
generally  not  feasible  unless  great  care  is  taken  to  match  the  distributions  in 
the  region  of  deep  fades.  As  a  result,  the  channel  must  be  carefully  modeled. 
Use  of  the  generalized  Rayleigh  distribution  has  encompassed  the  results  for  a 
wide  class  of  fading  distributions. 

Use  of  empirical  estimates  of  the  cross-correlation  moments  has  been  shown 
to  provide  highly  accurate  estimates  of  the  error  rate  unless  the  fading  is 
very  shallow.  In  addition,  it  has  been  shown  that  the  direct  use  of  moments 
based  on  integral  timing  offsets  between  the  codes  does  not  generally  provide 
good  estimates.  The  moments  must  be  calculated  for  arbitrary  timing  offsets; 
however,  the  moments  for  arbitrary  offset  can  be  readily  developed  from  the 
moments  for  integral  timing  offsets. 
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It  has  also  been  shown  that  the  effect  of  fading  on  the  error  rate  is 
significantly  reduced  for  highly  correlated  channels.  Further,  unless  the 
fading  correlation  is  either  close  to  unity  or  the  fading  on  the  intended 
signal  is  very  shallow,  long  spreading  codes  will  be  required  to  maintain  low 
error  rates. 


4.3.  Recommendations 

The  analysis  presented  here  has  been  restricted  to  the  case  of  slow, 
nonselective  fading.  With  this  basis,  extension  of  the  results  to  relatively 
fast,  nonselective  fading  would  provide  additional  insight  into  direct  sequence 
system  performance  over  generalized  Rayleigh  fading  channels.  Alternatively, 
the  effect  of  time  variations  in  the  channel  characterization  could  be 
explored,  that  is,  the  channel  fading  parameters  could  be  taken  as  random 
variables  with  some  assumed  distribution.  The  analysis  would  proceed  as 
presented  and  then  the  expectation  over  the  fading  parameter  distributions 
would  provide  the  long-term  average  error  rate.  In  general,  this  extended 
analysis  may  not  prove  tractable.  However,  the  single-term  estimates  of  the 
average  probability  of  error  could  be  analyzed  in  this  fashion  with  reasonable 
anticipation  of  success. 

An  alternate  area  of  investigation  would  be  the  analysis  of  diversity  gain 
(e.g.,  space  diversity)  possible  with  direct  sequence  systems  operating  over 
generalized  Rayleigh  channels.  This  analysis  would  be  particularly  important 
for  channels  whose  fading  depths  were  on  the  order  of  Rayleigh  fading  or 
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APPENDIX  A 

MATHEMATICAL  FUNCTIONS  AND  FORMULAS 

To  simplify  the  references  in  this  appendix,  they  are  abbreviated  to  two 
letters  designating  the  work,  followed  by  the  section,  page  or  equation  number 
in  the  reference.  The  two-letter  abbreviations  correspond  as  follows: 
AB-tAbramowitz,  1964],  ER-[Erdelyi,  1953],  GR-[Gradshteyn  and  Ryshik,  1965], 
HN-[Hansen,  1975],  MO-[Magnus,  Oberhettinger ,  and  Soni,  1966]  RA-[Rainville , 
I960],  VT-[Van  Trees,  1968]. 

1.  Pochhammer's  symbol  [AB6.1.22,  ER 1.21(5),  GRp.xliii,  HN2.4,  RAs.l8] 


(a)^  =  a(a+1 ) (a+2) . . . (a+k-1 )  =  T(a+k)  /  F(a)  ,  (a)Q  =  1  (A.1) 


2.  Hypergeometric  function  [ AB  15.1.1,  ER2. 1.1(2),  GR9.100,  HN10.9.2,  RAs .29 C 1 ) ] 


2F  (a,b;c;x)  =  I  ( c) 
k=0  '  ;k 


-  (a)k(b)k  xk 


(A. 2) 


3.  Linear  (Euler)  transformation  for  hypergeometric  function  [AB15.3.3,4; 
ER2. 1  .4(22), (23);  GR9.13K1);  HN10.9.2;  MOp.47;  RAs. 38(4) ,  (5)] 


2F1 (a,b;c;x)  =  (1  -  x)G-a-b  2F 1 ( c-a,c-b; c; x) 


(1  -  x)“  2F1(a,c-b;c;^‘) 


(A. 3) 
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4.  Hypergeometric  function  with  unit  argument  [AB15.1.20,  ER2. 1.3(14) 
HN7.4.14,  M0p.40 ,  RAs.32] 


2F1 (a,b; c; 1 )  =  [r(c)  r(c-a-b)]  /  [r(c-a)  r(c-b)] 
where  c  4  0 ,  -1 ,  -2 ,  . . .  ;  Re  (c-a-b)  >  0. 

5.  Confluent  (degenerate)  hypergeometric  function  [ AB 13-1.2, 

GR9. 210(1),  HN10.8.3,  MOp.262,  RAs.68(1)] 


i (a;b;x) 


00 


l 

k=0 


(a), 


(b)kk: 


k 


x 


6.  Kummer's  transformation  for  the  confluent  hypergeometric 
CAB  1 3 - 1  - 27 ,  ER6.3(7),  GR9. 212(1),  HN10.8.3,  MOp.267,  RAs.69(2)] 


1 F 1 (a;b;x) 


ex  1F1 (b-a; b; -x) 


7.  Generalized  hypergeometric  function  [ER4.1(1),  GR 9.14(1),  HN10.48.1 
RAs. 44 (2 ) ] 


F 

P  Q 


(a. 


.a- 


,ap’  brb2’  "  ' 


00 


l 

k=0 


(b,)k(b2)l...<bq)k  k! 


k 

x 


8.  Confluent  hypergeometric  function  of  three  variables  [HN48.21.18] 


*>2  ( a ,  b ;  c ;  x ,  y ,  z ) 


l  - 

k,n,m=0 


(a)k 

n!  m! 


(b) 

n _ 

(c), 

k+n+ra 


n  m 


,  GR9.122, 

(A. 4) 

ER6.K1), 

(A. 5) 

function 

(A. 6) 

,  MOp.62, 

(A. 7) 

(A. 8) 
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9.  Binomial  function  (binomial  theorem)  [ AB 15.1.8,  ER2.8(4),  HN1Q.4.6 
RAs.45(3)] 


(1  -  x)“a  =  ^FgCajjx)  =  2F1(a,b;b;x) 


10.  Binomial  coefficient  [AB6.1.21,  ER1.21(6),  M0p.4] 


x 

y  J 


=  r(x+i)  /  [ r(y+i )  r(x-y+i)]  , 


f  X  ]  =  (-1)n  (-x)  /  n! 

n  J  n 


11.  Generalized  Laguerre  function  [ ER 6.9.2(37) ,  MOp.336] 


(t) 


1F1 (-r;a+1;x) 


Lp(x) 


When  r  is  a  nonnegative  integer,  the  Laguerre  functions  are 
polynomials  (AB13.6.9,  ER6.9.2(36),  GR8. 972(1),  HN10.8.4,  RAs.112 

example , 

(x)  =  1,  L^a\x)  =  a+1-x,  L^a\x)  =  ^[x2-2x(a+2)  +  (a+1  )(a+2)] 

12.  Gamma  function  (Euler's  integral  of  the  second  kind)  [AB6.1;  ER1. 
MOp.1;  RAs.8,3.12] 


T(x+1)  =  /  dt  e-t  tX  ,  Re  x  >  0  ;  T(x+1)  =  x  T(x)  =  x! 

0 


,  MOp.33, 

(A. 9) 

(A. 10) 

(A. 11) 

Laguerre 
1)].  For 

(A. 12) 

;  GR8.31; 

(A. 13) 
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13.  Legendre's  duplication  formula  [AB6.1.18;  SRI. 2(15);  GR8. 335(1);  HN2.6,2.7; 
M0p.3,  RAs.19] 

F(2x)  =  22x_1  r(x)r(x+1/2)/r(1/2),  (2x)2k  =  22k  (x)k  (x+1/2)k  (A. 14) 

14.  Incomplete  gamma  function  [AB6.5,  ER9.1,  GR8.35,  HN6.10.1,  MOp.337] 

CO 

F(a,x)  =  /  dt  e-t  ta_1  =  7(a)  _  a-1  xa  ^ (a;a+1 ; -x)  (A. 15) 

x 

The  integral  representation  requires  larg  x I  <  w. 

15.  Beta  function  (Euler's  integral  of  the  first  kind)  [AE6.2,  ER1.5,  GR8.38, 
M0p.7,  RAs.16] 

B(a,b)  =  F(a)  F(b)  /  F(a+b)  =  /  dt  ta_1  (1-t)b_1 

0 

'n’/2  ~  .  ph  1 

=2  /  dt  (sin  t)  ~  (cos  t)^  ”  (A.16) 

0 

The  integral  representations  require  Re  a  >  0,  Re  b  >  0. 

16.  Incomplete  beta  function  [AB6.6,  ER2.5.3>  GR8.391,  HN 6.11 .12,  MOp.356] 

x 

B  (a,b)  =  /  dt  ta-1  (1-t)b_1  =  a"1  xa  _F  (a,1-b;a+1;x)  (A. 17) 

x  0  d  1 

17.  Normalized  incomplete  beta  function  [AB6.6,26.5;  ER2.5.3,  GR9.392] 
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I  (a,b)  =  B  (a,b)  /  B(a,b)  =  1  -  I  (b,a)  (A. 18) 

18.  Error  function  and  complementary  error  function  [AB7.1,  ER6.9.2(23), 

GR8.25,  HNp.507,  MOp.349,  RAp.127] 

x  2 

erf (x)  =  [2/r(1/2)]  /  dt  e_t  =  1  -  erfc(x)  =  2  P(21/2x)  -  1 

0 

=  2x  1F1 (1/2;3/2;-x2)  /  T( 1/2)  (A. 19) 

where  P  i3  the  standard  normal  cumulative  distribution  function. 

19.  Modified  Bessel  function  of  the  first  kind  (order  r)  [AB9.6.47, 
ER7.2.2( 12) ,  GR8.445,  HN10.7.11,  MOp.66,  RAs. 65(1)1 

Ir(2x)  =  xr  QF1 ( ; r+1 ; x2)  /  T(r+1)  (A. 20) 

20.  Hille-Hardy  formula  (bilinear  generating  function  for  Laguerre  polynomials) 
[ER10. 12(20),  GR8. 976(1  ),  HN48.21  .20,  MOp.242,  RAs. 120] 

l  tt+T)  Lka)(x)  Lka)(y)  =  (1-t)_1  r(a+1) 

k=0  la+i;k  K  K 

exp[t(x+y)/(t-1 )]  Ia[2(xyt)1/2/(1-t)l  ,  | tj  <  1  (A. 21) 


21.  Bilinear  generating  function  for  Laguerre  polynomials  [HN48.21.18] 
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i  tk  l<a><x)  L<  =  >(,)  =  ' 

k“0  (c+m-t-1)^  k  k+m  J  m! 


exp[t(x+y)/(  t-1 )]  $2t-m»c“a+n;c+1 ;y/(l-t) ,yt/( 1-t) ,xyt( 1-t)~2]  (A. 22) 


For  m  equal  to  zero  and  o  equal  to  a,  this  reduces  to  the  Hille-Hardy  formula. 

22.  Laguerre  series  for  modified  Bessel  functions  of  the  first  kind  [AB22.9.16, 
ER  10.12(18),  GR8. 975(3),  HN48.7.2,  MOp.242,  RAs. 113(2)] 


l  7-A-t-  (-Dk  t2k  Ja)(x2)  =  T(a+1)  (xt)'a  e-t  I  (2xt) 
k=0  U+a'k  *  a 


(A. 23) 


23.  Generalized  Marcum's  Q-function  [VTp.411] 


Qm(a,b)  =  /  dx  a  (x/a)m  exp(-  ^(x2  +  ac)]  I^tax) 


(A. 24) 


Q1  (a,b)  =  Q(a,b),  Qm(a,0)  =  1,  0^(0, b)  =  r(m,bV2)/r(m) 


(A. 25) 


24.  Marcum's  Q-function  [HN58.5.3,  VTp.395] 


Q(a,b)  =  expt-^(a2  +  b2)]  I  (a/b)k  I  (ab)  ,  alb 

k=0  K 


=  1  -  exp[-^(a2  +  b2)]  l  (b/a)^  I.  (ab) 

k=1  K 


b  1  a 


(A.  26) 
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25.  Integral  of  Bessel  function,  exponential  and  power  [AB11.4.28,  ER7.7.3(22), 
GR6.63K1)] 


7  p  ?  ,_i  (b/a)  r(V>  p 

J  dx  exp(-a  x  )  x3  I  (2bx)  =  -  .F.  t (r+s)/2;r+1 ; (b/a)  ] 

0  r  2a3  r(r+1)  1  1 

=  \  r[(s-r)/2]  (b/a)r  a-3  exp(b2/a2)  (~b2/a2)  (A. 27) 

p 

where  Re  a  >  0,  Re  (s+r)  >  0. 

26.  Interchange  of  order  of  summation  [HNs.4.1,  RAs.37] 


r* 

l 

n=0 


I  f(k,n) 
k=0 


oo 


l 

k=0 


00 


I 

n=0 


f (k,n) 


oo 


l 

n=0 


n 

Y  f(k,n-k) 
k=0 


oo 


l 

n=0 


n 

l  f(n-k,k) 
k=0 


27.  Powers  of  series  [AB3.6.20,  GR0.314] 


(A. 28) 


(  l  ak  x*  )“  =  I  ck  x  ’ 

k=0  *  k=0 


(A. 29) 


where 


c„  =  a”  ,  c_  =  TT"  l  (kn-m+k)  a k  cm_k  for  m^l 


'0  0  ’  m  "  ma 


(A. 30) 


0  k=1 


28.  Double  factorial  function  [ABp.258,  GRp.xliii] 
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(2n) ! !  =  (2)(4)...(2n)  =  2n  n!  ,  (2n-1 ) ! !  =  (1 )  (3) . . • (2n-1 )  (A. 3D 

29.  Generating  function  for  Stirling  numbers  of  the  first  kind  [AB24.1.3, 

HN52.2.1] 

(-x)  =  (-1  )n  l  xm  S^n.m)  (A. 32) 

m=0 

30.  Recurrence  relation  for  Stirling  numbers  of  the  first  kind  [AB24.1.3] 

S1  (n+1  ,m)  =  S 1  ( n , m- D  -  nS^n.m)  ,  n>m21  ;  S^n.O)  =  «0n,  (A. 33) 

S^(n,D  =  (-Dn  1  (n— 1 ) ! ,  S^n.n)  =  1,  S^(n,n-1)  =  -n(n-1)/2  (A. 34) 

31.  Generating  function  for  Stirling  numbers  of  the  second  kind  [AB24.1.4, 

HN52.2.15] 

(_x)n  =  l  (-1)m  <x)B  S2(n,a)  (A. 35) 

m=0 

32.  Recurrence  relation  for  Stirling  numbers  of  the  second  kind  [AB24.1.4] 

S2(n+1,m)  =  mS2(n,m)  +  S2(n,m-1)  ,  n2m2.1  5  (A. 36) 

S2(n,0)  =  6^  S2(n,1)  =  S2(n,n)  =  1,  S2(n,n-1)  =  n(n-1)/2 


(A. 37) 
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33.  Series  involving  Stirling  numbers  of  the  second  kind  [HW6.7.22] 

n  m 

7  (-D  (-n).  km  /  k!  =  2n  y  (-n).  (-2)"*  S.(m,k)  (A. 38) 

k=1  k=1 

34.  Generating  function  for  generalized  Euler  numbers  of  order  (-m)  [HN51.7, 
MOp.32] 


(cos  t)m  =  Y  (-1)kt2k  E»~m)  /  (2k)! 

k=0 


(A. 39) 


From  Hansen  [HN88.1.2],  we  have 


(cos  t)  =  l  (2k)I 
k=0  K  *’■ 


t2k  [2-m  “  (■  » 

n=0  *• 


(2n-m)2k] 


(A. 40) 


By  comparing  Eq.(A.40)  and  Eq.(A.39),  we  obtain  [Hanlon,  Peterson,  and  Gardner, 
1980] 


(-m)  _  -m  ? 

2k  ‘  Ln 

n=0 


r  \ 

m 

n 


(2n-m) 


2k 


(A. 41) 


35.  Generating  function  for  Bernoulli  numbers  [AB23.1.1,  ER 1.13(1),  GR9.610, 
HN50.5.3,  RAs. 153(8)] 


oo 

t  /  (efc  -  1)  =  l  Bn  tn  /  n!  ,  | t| <2W 


(A. 42) 


n=0 


For  example, 
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Bq  =  1  ,  B1  *  -1/2  ,  B2  =  1/6  ,  B4  =  -1/30 


(A. 43) 


36.  Series  expansion  for  logtcos  x)  in  Bernoulli  numbers  [AB4.3.72,  ER1.20(S), 

GR1. 518(2),  HN50.6.8] 


log( cos  x)  =  l  (-1)n  22n-1  (22n-1)  B2n  x2n/[n(2n) !  ] ,  lxl<ir/2  (A. 44) 
n=  1 
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APPENDIX  B 

PROBABILITY  DENSITY  FOR  FADING  INTENSITY 

B.l.  Univariate  Density 

Let  y1 , . . . ,y  be  independent  and  identically  distributed  normal  random 

2 

variables  with  mean  P  and  variance  o  .  The  Euclidean  norm  of  these  variables, 
designated  f,  is  defined  by 

f  =  (  l  y\  )1/2  .  (B.D 

i=1  1 

It  is  well-known  that  f  has  a  noncentral  chi  distribution  with  n  degrees  of 
freedom  and  its  probability  density  is  given  by  [Miller,  Bernstein,  and 
Blumenson,  1958] , [Nakagami ,  I960] 

_  fl!  _  . 

Pf(x)  =  — exp!-m(x  +  B'/n)  Im_i  (2m6x/fi)  ,  x  >  0  (B.2) 

ns™- 

where  I  is  a  modified  Bessel  function  of  the  first  kind  of  order  r  and 
r 

o  *  /  p 

m  =  n/2  ,  n  =  no,  8  =  n  u  •  (B.3) 

Although  this  development  requires  n  to  be  a  positive  integer,  Eq.(B.2)  is  a 
valid  fading  probability  density  for  all  m  not  less  than  1/2.  In  this  context, 
f  is  referred  to  as  a  generalized  Rayleigh  variate  (generalized  n-distributed 
variate)  [Nakagami,  I960] .[Miller,  Bernstein,  and  Blumenson,  1958].  For  a 
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random  variable  whose  probability  density  is  given  by  Eq.(B.2),  its 
distribution  is  represented  by  f  ,vR(m,6,Q).  The  moments  of  f  are  readily 
evaluated  by  using  Eq.(A.27)  and  are  given  by  [Nakagami,  I960] ,[Krishnan,  1967] 


V  V /?  ? 

<f  >  =  r(m  +  v/2)  (fi/m)  1F1 (-v/2;m; -mBVft)  /  r(m) 

=  r( v/2  +  1)  («/m)V/2  L^"1  ^-me2/^)  ,  v+2m  >  -1  (B.4) 

where  <  >  represents  expectation,  ^  is  a  confluent  hypergeometric  function, 
(  S  ) 

and  L  is  a  generalized  Laguerre  function.  Note  that  the  mean  square  value 
r 

of  f  is  just 

<f2>  =  0  +  62  .  (B.5) 


The  generalized  Rayleigh  distribution  encompasses  a  wide  class  of 

distributions  which  are  useful  in  the  study  of  communications  over  fading 

channels.  For  m  equal  to  one  and  0  equal  to  zero,  f  has  a  Rayleigh 

distribution.  For  m  equal  to  one  and  8  arbitrary,  f  has  a  Rician  distribution 

2 

(Nakagami-Rice  or  n-distribution) .  In  this  case,  6  is  the  specular  power  and 
^  is  the  mean  square  value  of  the  Rayleigh  component.  The  depth  of  fading  can 
be  measured  by  the  specular  power  ratio  (the  ratio  of  the  specular  power  to  the 
mean  square  value) ,  designated  R  and  defined  by 
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R  =  g2  /  ( n  +  B2)  .  (B.6) 

The  depth  of  fading  increases  with  decreasing  R.  However,  Rician  fading  can 

never  be  deeper  than  Rayleigh  fading  (S  =  0).  For  3  equal  to  zero  and  m  not 

less  than  1/2,  f  has  a  Nakagami  distribution  (m-distributed) .  In  this  case, 

2  2 

the  mean  square  value  of  f  is  J!  and  the  variance  of  f  is  ( Q  /m) .  Hence,  m  is 

2 

the  inverse  of  the  normalized  variance  of  f  [Nakagami,  I960].  The  depth  of 
fading  increases  with  decreasing  m.  For  values  of  m  less  than  one,  the  fading 
is  more  severe  than  Rayleigh  fading,  and  for  values  of  m  greater  than  one,  is 
shallower  than  Rayleigh  fading.  Consequently,  Rayleigh  fading  i3  a  special 
case  of  both  Rician  and  Nakagami  fading.  In  the  general  case,  the  fading 
decreases  as  either  m  or  6  increases. 

The  cumulative  distribution  function  for  a  generalized  Rayleigh  variate  is 
given  by 


F„(x)  =  1  -  Q[(2m32/.Q)1/2  ,  (2mx2/S2) 1/2]  (B.7) 

i  m 

where  Q  i3  a  generalized  Marcum's  Q- function.  To  obtain  an  alternate 

m 

expression,  we  make  use  of  the  Laguerre  series  for  modified  Bessel  functions 
(Eq.(A.23))  to  express  the  probability  density  as 

00  /  1  \  ^ 

p~(x)  =  exp(-m3“/fl)  7  — r~ 

f  nr  On)  k=o  (m)k 

Integration  of  the  probability  density  is  then  straightforward  and  provides 


mx 


m+k-l 


.  ( m- 1 ) ,  2  i  \ 

Lk  (mg  /$:) 


C  B .  S ) 
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[Hanlon  and  Gardner,  1980] 


Ff(x)  = 


r(m+1) 


exp(-mp  /ft)  l 
k=0 


(-1) 
(m+1 ). 


2 

mx 

Q  j 


m+k 


(B.  9 ) 


The  probability  that  one  generalized  Rayleigh  variate  exceeds  another  when 
they  are  independent  is  given  by 


Pr{fp  >  f  }  =  /  dx  p_(x)  /  dy  pf(y) 

0  2  0  1 


/  dx  p  (x)  F_(x) 

0  r2  1 


<F  (f  )> 
r1  * 


(B. 10) 


By  using  Eq.(B.9)  and  Eq.(B.4),  we  readily  obtain 


Pr{f2  >  s  exp(-n18^/ni)  l  (-1)k  (ny'ym^)  1 


m,  +k 


(01,-1)  2  (m_-1)  2 

Lk  (m1 61  /n!)  Lk+m1  (-®2B2/?'2)  • 


(3.11) 


Note  that  if  the  ratio  (m^^/m^G^)  is  greater  than  one,  the  series  does  not 
converge  and  the  representation  given  by  Eq.(B.II)  is  formal.  In  this  case,  we 
can  use  the  alternate  expression 


Pr{f2  >  f^  =  1  -  Pr { f  1  2  f2J 


(B. 12) 
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which  can  be  written  from  Eq.(B.II)  by  inspection.  Several  special  cases  of 
the  general  result  are  of  interest: 

(a.)  For  f i  and  f^  both  Rician  (m^m^l),  use  of  Eq.(A.22),  Eq.(A.20),  and 
Eq.(A.26)  provides  [Schwartz,  Bennett,  and  Stein,  1966],  [Van  Trees,  1968] 


(HUSg/cnj+n  )]‘ 


[2B^/(:i1+22)]1/2) 


[2.1/(ai+Sl2)]  exp[-(62+S2)/(ni+a2)]  IQ[2S1S2/(ni- 


(3.13) 


(b.)  For  f1  and  f2  both  Nakagami  m-distributed  ( 3 1  =  3^=0 ) ,  Eq.(B.II) 
readily  reduces  to 

r  p  )  m, 
m1  “5  1 

2F1(m1  ,m1+m2;m1  +  1;-m122/m221)  .  (B.14) 

After  linear  transformation  of  the  hypergeometric  function  and  use  of  Eq . ( A . 1 7 ) 
and  Eq . ( A . 1 8 ) ,  we  obtain 


Pr{f2  >  f^  = 


m^  +m  -1 


It(nrm2) 


t  =  [m. 


/  (m.2„ 


(3.15' 


(c.)  For  f2  one-sided  Gaussian  (m2  =  1/2,  S 2  =  0),  we  obtain  ty 


subst itution , 
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Pr{f2  >  f^ 


m.,-1/2 

“l 


“  (m  +1/2), 

exp(-»,ef^,)  jo~ 


2m 


m^+k 


ht  • 


(B.16) 


Table  B.1  lists  the  results  for  several  additional  special  cases. 


B.2.  Bivariate  Density 

Let  the  pairs  (y^  1  ,y12)  *  •  ♦  • * (yni  be  indePendent  and  identically 

distributed  where  the  elements  and  y^2  are  jointly  normal  for  all  j  with 

2  2 

means  u1  and  u2,  variances  and  and  correlation  coefficient  P.  Then  the 
norms,  f1  and  f2,  defined  by 

ft  ■  (  I  y2u  )1/2  (b.  17) 

1  J=1 

each  have  a  noncentral  chi  distribution  with  n  degrees  of  freedom  and  the 
bivariate  density,  for  the  magnitude  of  P  less  than  one,  is  given  by  [Miller, 
Bernstein,  and  Blumenson,  1958] ,[Krishnan ,  1967] 


1  -m  2  2 

pf  f  (x1tx2)  =  A  (pD.|D2)  x1  x2  expC-B^-B^) 


k=0 


Ck  Ir(2D1x1)  Ir<2D2x2)  yZpBj'V/2*^) 


(B. 18) 


where 
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TABLE  B.1 

SPECIAL  CASES  FOR  PROBABILITY  THAT  ONE  GENERALIZED  RAYLEIGH  VARIATE  EXCEEDS 
ANOTHER  WHEN  THEY  ARE  INDEPENDENT 


f1 

f2 

(m1,81;m2,S2) 

Pr{f2  > 

Rician 

Rayleigh 

(1, 8^1,0) 

Cn2/(^1  +  n2)]  expt-B2/^  +  «2)] 

Nakagami 

Nakagami 

(m1 ,0;m2,0) 

It(m1,m2)  ,  t  =  [m1n2/(m1n2+m2n1)] 

Nakagami 

Rayleigh 

(m^ ,0; 1 ,0) 

ml 

[m^y^-m^)]  1 

Nakagami  'Gaussian 

(a^.0;  1/2,0) 

It(n1(1/2)  ,  t  •  [2m1^2/(«1+2m1«2 

Rayleigh 

Rician 

d,0;i,62) 

i  -  [<y(n1+n2)]  exP[-B2/ayo2)] 

Rayleigh 

Nakagami 

O 

tf* 

o 

nu 

1  -  [m^^+m^)]  41 

Rayleigh 

Rayleigh 

( 1 ,0; 1 ,0) 

[sycsyiy] 

Rayleigh  *Gaussian 

( 1 ,0; 1 /2 ,0) 

i  -  [r1/(o1+2o2)]1/2 

•Gaussian 

Nakagami 

(1/2,0;m2,0) 

It(1/2,m2)  ,  t  =  Cn2/(n2+2m2fi1)] 

•Gaussian 

Rayleigh 

(1/2,0; 1 ,0) 

[S22/(n2+201  )31/2 

•Gaussian  *Gaussian 

(1/2,0; 1/2,0) 

f  sin"1[02/(01+02)]1/2 

•Gaussian  refers  to  one-sided  Gaussian. 
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(2m-2).(m). 

'  ck  -wfiT  <B-’9> 

1/21/2  1/21/2 

ti>1  -  8^^  -  p82^-)  ^2  »  ^2  =  “  P^i^i  (B.20) 

Bt  =  mtfljO-p2)]"1  ,  Di  =  m^t^^d-p2)]*1  (B.21) 

A  =  4B1B2(1-p2)r(m)  (mfi^1/2n21/2)ni~1  expMmB*/^ )  -  (D^/B^]  .  (B.22) 

Note  that  in  the  Rician  case  (ms  1),  we  have 


lim  C.  s  lia  (2m-2),  /  (m-1).  s  g  (B.23) 
rn-,-1  *  m-t-1  K  K  K 

where  is  Neumann's  factor,  which  is  defined  to  be  equal  to  2  when  k  is  not 
zero,  and  equal  to  1  when  k  is  zero.  As  before,  Eq.(B.l8)  is  a  valid  bivariate 
probability  density  for  all  m  not  less  than  1/2  and,  in  this  context,  is  the 
joint  probability  density  for  correlated  generalized  Rayleigh  variates  with 
marginal  distributions  f1  ^  R(m, B^S^)  and  f2  ^  R(m,S2,fi2). 

In  general,  the  expression  for  either  the  joint  moments  or  the  correlation 
between  f1  and  f2  is  very  complicated  [Krishnan,  19673;  however,  for  the 
Nakagami  or  Rayleigh  cases,  the  normalized  covariance  of  f2  and  f^  (power 
correlation  coefficient)  is  very  simply  related  to  the  correlation  coefficient 
of  the  underlying  Qaussians,  that  is,  [Nakagami,  I960] 


V 


Cov(f^,f2)  /  [Var(f^)  Var(f^)]1/2  =  p2  . 
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(B.24) 


The  probability  that  one  generalized  Rayleigh  variate  exceeds  another  in 
the  correlated  case  is  given  by 

x. 

oo  2 

Pr{f_  >  f.)  =  /  dx_  /  dx.  p-  (x1tx5)  .  (B.25) 

*  0*0  1r2  '  * 

The  inner  integration  can  be  evaluated  by  expanding  the  exponential  and  Bessel 
functions  involving  x1  (Eq.(B.l8))  into  nested  power  series  in  the  variable 
x1  and  interchanging  the  order  of  integration  and  summations.  The  outer 
integration  is  then  readily  evaluated  by  Interchanging  the  order  of  integration 
and  summations  and  using  Eq.(A.27).  Use  of  Eq.(A.28)  and  Eq.(A.II)  then  yields 


Pr{f  >  f  }  =  (1-p2)m  expC-mr/n.)  \  (SI  /n  )k+m  (-1)1 

*  k=G  * 


Gk(m,D2/BrD2/B2,P2) 


(B.26) 


where 


k  (2m-2)  k-n  (1-k-m)  (1+k+m) 

Gk(m,x,y,t)  *  J  (a-1)  \\  (txy)"  \  s!  * 


a  ,  (nwa-1 ) ,  .  (n-wn-1 )  ,  , 

1  tc-n-s  Lk+nM>s  {  y)  • 


(B.27) 
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As  in  the  independent  case,  if  the  ratio  (Gg/Oj)  is  greater  than  one,  we  can 
evaluate  the  complement  of  the  probability  (Eq.(B.12)). 

For  the  case  where  f^  and  f2  are  independent  (psO),  the  only  contribution 
from  the  sums  over  s  and  n  in  is  the  first  term  (ns0,s=0)  and  Eq.(B.26) 
reduces  to  Eq.(B.II)  for  m^  and  m2  equal  to  m.  The  result  also  simplifies 
considerably  if  f^  and  f2  are  both  Nakagami  m-distributed.  In  this  case, 
and  D2  are  both  zero  and  the  only  contribution  from  the  sum  over  n  in  is 
the  first  term  (n=0).  By  noting  that 


we  then  have 


(B.28) 


Pr{f2  >  f^ 


fz-ii  ,,  2,»  ;  n  (">k  .  ,.k ..  ..  , 

l»  J  >  Jz0  kf  («!),;  'W 


k+m 


2F1(-k,k+2m;m;p  )  . 


(B.29) 


If  we  further  specialize  to  the  case  where  f^  and  f2  are  both  Rayleigh 
distributed  (m*1),  we  obtain 


Pr{f2  >  f1 >  = 


(1-p2) 


i, 


(-I)1'4’1  gF^I-k.k+ljljp2) 


(B.30) 
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APPENDIX  C 

SERIES  EXPANSION  FOR  PROBABILITY  OF  ERROR 


The  probability  of  error  ia  given  by 


®  w 

Pr{E}  =  /  dy  /  dx  p  (x,y)  .  (C.1) 

0  0  sz 

To  develop  a  series  expansion,  the  joint  density  of  S  and  Z  is  expanded  in  a 
one-dimensional  Taylor  series  about  x  equal  to  a,  that  is, 


=  I  **=■>. 

3Z  n=0  n! 


— n  psz(x,y) 

3xn  sz  x=a 


(C.2) 


Hence,  the  probability  of  error  can  be  represented  as 


Pr,E1  ■  o'  dy  I  *  Jo  W 


(C.3) 


By  interchanging  the  order  of  integrations  and  summation  and  by  noting  that  the 
partial  derivatives  evaluated  at  x  equal  to  a  are  not  a  function  of  x,  we 
obtain 


Pr{E}  * 


00 


x*a 


/  dx  (x-a)n  . 
0 


(C.4) 


After  evaluating  the  integration  over  x,  we  have 


Pr{E}  =  l  7-^rrr  /  dy  p  (x,y)  I  [(y-a)n+1  -  (-a)n+1],  (C.5) 

n=0  0  3xn  sz  x=a 


.  n+1  • 


By  using  the  binomial  theorem  (Eq.(A.9))  to  expand  (y-a)n+1 ,  Eq.(C.5)  becomes 

®  00  JJ  |  11^1  t  \ 

PrW  ■  J„  rzm  I «  J,  f”"1) y"  •  <c-6> 

After  again  interchanging  the  order  of  integration  and  summation  and  then 
interchanging  the  order  of  integration  and  partial  differentiation,  we  obtain 


?rW  *  1  ^  1  (”«l)  S  l  ^  '  <C- 


7) 


The  joint  density  of  S  and  Z  equals  the  product  of  the  conditional  density  of  Z 
given  S  times  the  marginal  density  of  S;  Eq.(C.T)  then  yields 


*  Jo  ^7  f  m 


~  P,(x)  /  dy  ym  P,|.(y|s=x)  |  . 


3  X 


Z  3 


(C.8) 


x=a 


Since  the  probability  density  of  Z  is  symmetric  about  the  origin,  the  remaining 
integral  is  half  the  m-th  conditional  absolute  moment  of  Z  given  S,  that  is, 


< ! Z |*  S«x>  »  2 


/  dy  y®  p  |  (y|S=x) 
0  z's 


(C.9) 


where  <  >  represents  expectation.  After  translating  the  index  of  summation,  we 
obtain  the  series  expansion  for  Pr{E}  about  x  equal  to  a: 


9M 


■  j  ni,  4r  j,  (») 


.  n-1 


3  x 


n-1 


Ps(x)  <|z|m 


S=x> 


(C.10) 


x»a 


For  the  case  where  a  equals  zero  (Maclaurin  series) ,  the  only  contribution  from 
the  inner  series  is  for  m  equal  to  n.  In  this  case,  we  have  [Hanlon  and 
Gardner,  1979  and  1980] 


Pr(E)  = 


.  tn-1 


3  x 


i-l 


(C.11) 


For  the  special  case  where  S  and  Z  are  independent,  the  conditional  absolute 
moments  of  Z  are  just  the  absolute  moments  of  Z  and  Eq.(C.IO)  becomes 


««  ■  i  l  TT  !,  ' 


n  (-a)n-m  <  |z  r> . 

m 


(C. 12) 


Equation  (C.11)  becomes 


pr[E}  =  j  J  Jj-  p^n_1)(0)  <|z|n>  .  (C.  13) 

ns  1 

It  should  be  noted  that  in  Eq.(C.12)  the  average  power  of  the  interference, 
i.e.,  <Z2>,  contributes  to  the  second  and  higher  terms  of  the  series  over  n. 
In  Eq. (C. 13) ,  the  average  power  of  the  interference  contributes  only  to  the 
second  term  of  the  series.  Since  one  would  intuitively  expect  the  probability 
of  error  to  be  strongly  dependent  on  the  average  power  of  the  interference, 
Eq.(C.13)  would  appear  to  be  a  more  natural  ordering  likely  to  provide  faster 
convergence.  Additionally,  the  derivatives  of  the  intended  signal's 


A 
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probability  density,  p^n_1\  is  in  general  of  simpler  form  when  evaluated  for  a 

s 

equal  to  zero,  as  opposed  to  some  other  value  of  a.  Similarly,  Eq.(C.II)  is 
generally  preferable  to  Eq.(C.IO). 


As  an  example,  the  series  expansion  given  by  Eq.(C.13)  is  used  to  evaluate 
the  probability  of  error  when  S  and  Z  are  independent,  S  is  a  generalized 
Rayleigh  variate  with  distribution  parameters  (n^ ,  g^,  where  2m1  is 
restricted  to  be  an  integer,  and  |z  |  is  a  generalized  Rayleigh  variate  with 
parameters  (m,,  ^ *  From  Eq.(B.8) ,  it  i3  readily  shown  that  the  even 
derivatives  of  the  density  of  S  at  the  origin  are  zero  and  the  odd  derivatives 
at  the  origin  are  nonzero  only  if  the  order  of  the  derivatives  is  greater  than 
or  equal  to  (Zm^l).  The  derivatives  are  then  (for  nlm^) 


(m.  -1) 


(On— A  ^  _  o  v*“i  9 

Pg  1;(0)  =  (-1)  [(2n)  !/n!  ]  (m,/^)"  expC-m^*/^)  (m^/flj)  (c.14) 

The  moments  of  |Z|  are  given  by  Eq.(B.4)  and  the  probability  of  error  is 
readily  shown  to  be  [Hanlon  and  Gardner,  1980] 


i  ?  n+m. 

Pr{E}  =  j  exp(-m18^/01 )  £  (-1)n  (m^/m^) 


(m.-l)  „  (m.-l)  9 

Ln  (nlW  Vm,  • 


This  result  is  consistent  with  Eq.(B.II). 


(C. 15) 
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APPENDIX  D 

CROSS-CORRELATION  MOMENTS  FOR  RANDOM  CODES 


D.l.  Momenta  Ion  Integral  Offset 

The  probability  mass  function  for  the  code  cross-correlation  under  the 
condition  of  integral  offset  for  binary  symmetric  random  codes  of  length  L  (L 
chips  per  data  bit)  is  given  by 


Pr{c  s  k} 


-L  1  k  1  L 
k  mod2  =  L  mod2 


0 


,  otherwise. 


(D.l) 


It  can  be  readily  shown  that  the  odd  moments  are  zero  and  the  even  moments  are 
given  by 


<c2n>  =  2"L  \  |j|j  (2k  -  L)2"  -  E^L)  (D.2) 

where  Eg”^  is  a  generalized  Euler  number  of  order  (-L).  However,  this  form  is 
not  very  convenient  since  L  can  be  very  large.  By  using  the  binomial  theorem 
(Eq.(A.9)),  Eq.(A.38),  and  the  generating  function  for  Stirling  numbers  of  the 
first  kind  (Eq.(A.32)),  the  even  moments  can  be  expressed  as  a  polynomial  of 
degree  n  in  the  code  length,  i.e.. 
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<c2n>  =  J  A  Lm  for  n  >  1  (D.3) 

m=1  m’n 

where  the  coefficients  (A  )  are  given  in  terms  of  Stirling  numbers  of  the 

m,n 

first  (S^)  and  the  second  (S^)  kinds  by  the  relation  [Hanlon,  Peterson,  and 
Gardner,  1980] 


A 


m,n 


2n-m  m-1 

7  2k  T  (-1;  S. (2n-r-k,m-r)  S.(2n-r ,2n-r-k) 

k»0  r»0  1  ' 


(D.4) 


Provided  n  is  not  too  large,  this  expression  readily  generates  the  polynomial 
coefficients.  However,  it  is  desirable  to  obtain  closed- form  expressions  for 
the  coefficients  of  at  least  the  higher-order  powers  of  L,  3ince  for  large 
values  of  L  (relative  to  n  ) ,  the  moments  can  be  approximated  very  accurately 
by  one  or  two  higher-order  terms  of  the  polynomial.  To  obtain  the  closed-form 
expressions  for  the  coefficients,  we  make  use  of  the  characteristic  function 
for  the  cross-correlation.  By  using  the  probability  mass  function,  we  have 


<eJ“c> 


0j<u(2k-L) 


(cos  to) 


L 


(D.5) 


By  applying  the  moment  theorem  and  Taylor's  theorem  to  Eq.(D.5),  we  obtain  for 
the  coefficients: 


2n 


A  =  (-1  )n  -1=-  ■  d-*  -[log (cos  u)]m 
m,n  ml  d(jJ2n 


0>sO 


(D.6) 


By  applying  the  series  expansion  for  log(cos  ot)  (Eq.(A.44))  and  the  result  for 
powers  of  series  (Eq.(A.29)),  we  have 
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A 


m,n 


(-1) 


m+n 


(2n)l  H 
m!  n-m,m 


where 


0,m 


=  2 


-m 


k,m 


_2_ 

k 


k 

l 

r=1 


(rm-k+r)  G  H. 


r  k-r,m 


(D.7) 


(D.8) 


G  =  (-1)r  22r+1  (22r+2-1)  B_  5  /  t (r+1 ) (2r+2) ! ]  ,  (D.9) 

r  2r+2 

and  the  Bn  are  Bernoulli  numbers.  For  example,  for  the  highest-order 
coefficients  we  have 


A  =  (2n) !  /  (n!  2n)  =  (2n-1)!! 

n,n 


(D. 10) 


An-1  ,n 


-n(n-1 )  A  /  3 
n  ,n 


(D.11) 


An-2 ,  n 


-(n-2)(5n+1)  An  /  30  . 
n-l  ,n 


(D . 12 ) 


2 

For  large  values  of  L  relative  to  n  ,  the  even  moments  are  approximately  given 
by 


<c2n>  =  (2n-1 ) ! I  Ln  [1  -  n(n-1)/3L]  . 


(D. 13) 


This  expression  is  exact  for  the  second  and  fourth  moments.  Table  D.1  lists 
the  coefficients  for  the  first  five  even  moments. 
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TABLE  D.  1 

COEFFICIENTS  FOR  POLYNOMIAL  EXPRESSIONS  OF  CROSS-CORRELATION  MOMENTS  FOR  RANDOM 


CODES  UNDER 

INTEGRAL 

OFFSET,  <c2n> 

n 

=  l  A  n 

m=1  m'n 

n 

A„  „ 
n,n 

An-1,n 

An-2 , n  An-3 , n 

An-4,n 

1 

1 

2 

3 

-2 

3 

15 

-30 

16 

4 

105 

-420 

588  -272 

5 

945 

-6300 

16380  -18960 

7936 

D.2.  Momenta  for  Arbitrary  Offset 

The  moments  for  an  integral  offset  can  be  used  to  generate  the  moments  for 
an  arbitrary  offset.  From  Eq.(2.5),  by  talcing  the  expectation  over  x,  which  is 
uniform  on  the  interval  (0,1),  while  holding  cQ  and  c1  fixed,  we  obtain 


..n 

<c 


W  = 


n+1 


n 

L 

m=0 


m  n-a 

=0  °1 


(D.14) 


Since  Cq  and  are  independent  for  random  codes  and  the  odd  moments  of  c  are 
zero,  then  the  odd  moments  of  c  are  zero  and  the  even  moments  are  given  by 


(D.  15) 


<c2n>  = 


2n+1 


l  <c2m><c2n“2m> 


_ 1_ 

s  2n+1 


n 

I 

m=0 


?(-L) 

"2m 


E 


(-L) 
2n-2m  * 


It  is  clear  from  Eq.(D.3)  and  Eq.(D.15)  that  the  even  moments  of  c  are  also 
polynomials  of  degree  n  in  the  code  length  and  can  be  expressed  as 

n  _ 

<c  >  =  y  Lm  for  n  >  1  .  (D.16) 

m=1  m’n 

By  solving  for  the  coefficients  in  terms  of  the  A  ,  we  obtain  [Hanlon, 

m,n 

Peterson,  and  Gardner,  1980] 


min(m,k) 

I 


m’n  2n+1  k=0  r=max(0,k-n+m)  r’k  m-r’n-k 


(D.17) 


where  we  define  ^  =  6^^,  that  is,  the  Kronecker  delta.  For  large  values  of 
2 

L  relative  to  n  ,  we  can  approximate  the  even  moments  by  two  terms  of  the 
polynomial  involving  the  highest-order  powers  of  L,  that  is, 


<s2n>  =  -fnTT1--  Jo  (;)  (S)’1  I'  -  [»2*(n-m)2-n]/3U  .  (0.18) 

This  expression  is  exact  for  the  second  and  the  fourth  moments.  Table  D.2 
lists  the  coefficients  for  the  first  five  even  moments.  It  should  be  noted 
that  the  moments  for  arbitrary  offsets  differ  significantly  from  those  for 
integral  offsets. 


TABLE  D.2 


COEFFICIENTS  FOR  POLYNOMIAL  EXPRESSIONS  OF  CROSS-CORRELATION  MOMENTS  FOR  RANDOM 


CODES, 

<c2n>  =  l 

Dn  n  L“- 

ms  1 

in  f  n 

n 

D 

D  , 

D  „ 

D 

D 

n,n 

n-1  ,n 

n-2,n 

n-3,n 

n-4,n 

1 

2/3 

2 

7/5 

-4/5 

3 

36/7 

-64/7 

32/7 

4 

249/9 

-912/9 

1212/9 

-544/9 

5 

2190/11 

-13680/11 

34152/11 

-38528/11 

15872/11 

102 
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